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Preface 



This dissertation intends to present a unified view of my researcli in tlie area 
of integrable systems, which has been carried out over the last four years 
under the insightful supervision of Profs. F. Finkel, A. Gonzalez-Lopez and 
M. A. Rodriguez. Particular emphasis has been laid on the connection be- 
tween (partially) solvable spin models of Calogero-Sutherland (CS) type and 
Haldane-Shastry (HS) chains, which can be rightly regarded as the leitmotif 
of this thesis. 

The dissertation is divided into five chapters. In the first one we sum- 
marize some aspects of the theory of CS and HS models that are needed in 
later chapters and introduce some notation. Related original results are in- 
troduced when appropriate. In Chapter [2] we rigorously analyze the CS/HS 
connection, which enables us to compute the partition function of the BCn 
HS chain in closed form. In Chapter [3] we turn our attention to quasi-exactly 
solvable (QES) generalizations of the standard CS models presenting only 
near-neighbors interactions. A systematic explanation of their quasi-exact 
solvability is given in terms of known invariant flags. In Chapter H] we make 
use of the results in the two preceding chapters to analyze the algebraic 
spectrum of a novel spin chain with nearest-neighbors interactions which is 
naturally related to a QES Schrodinger operator of CS type. This work ends 
with a short summary and a concluding note on future projects. 

The most remarkable contributions presented in this thesis are the fol- 
lowing: 

(i) We complete the first mathematically rigorous treatment of the CS/HS 
connection. Our approach hinges on several convergence results that 
make Polychronakos's "freezing trick" precise. 

(ii) We provide the first systematic analysis of spin QES models of CS 
type presenting near-neighbors interactions, thereby obtaining several 
families of totally explicit eigenf unctions. 

(iii) We define an HS chain with nearest-neighbors interactions that presents 
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some of the features of a QES system. The study of its algebraic states 
rehes on a nontrivial refinement of the freezing trick and on the exphcit 
resuhs obtained in Chapter [31 

(iv) We introduce a natural Hamiltonian system admitting the maximum 
number of independent first integrals, just as the original Calogero and 
hyperbolic Sutherland models. This is the first example of a maximally 
superintegrable system in an A^-dimensional space of nonconstant cur- 
vature. 

A significant part of the material presented in this memoir has been 
adapted from Refs. [15,63-70]. The presentation of the main results in Chap- 
ters [2] and m however, greatly differs from that of the original articles. In 
particular, most of the rigorous proofs provided in these chapters appear in 
printed form here for the first time. 



Chapter 1 

Calogero— Sutherland models 



1.1 Introduction 

The discovery of Calogero-Sutherland (CS) models [35, 186, 187] is a land- 
mark in the theory of integrable systems which has exerted a far-reaching 
influence on many different areas of Mathematics and Physics. The intrinsic 
mathematical interest of these models and their classical analogs [153, 157] 
lies in their manifold connections with such disparate topics as representation 
theory [22,72,80,177], the theory of special functions and orthogonal poly- 
nomials [14,60,190,191], harmonic analysis [40,108,174], and complex and 
algebraic geometry [71,114,119,154,202]. On the other hand, CS models 
have found a wide variety of applications to diverse areas of Physics, in- 
cluding quantum field and string theory [8,54,85,95,149,151], the quantum 
Hall effect [11,127,169], soliton theory [1,33,128,133], fractionary statis- 
tics [37,99,163] and random matrix theory [38,181]. 

In this chapter we summarize the basic properties of CS models that 
shall be needed in the sequel. We prefer not to follow a historical approach 
to this subject, but rather to develop a modern bird's-eye presentation of 
the essentials of the theory. Thus the most advanced machinery becomes 
available from the very beginning, while the interest of the reader is not 
suffocated by unnecessary details. Our survey of CS models does not intend 
to be exhaustive. In particular, we have omitted several aspects of the theory 
(e.g., spin CM models [121,140,142,204], Ruijsenaars-Schneider systems and 
g-Dunkl operators [9,13,43,141,176,198], supersymmetry [48,51,93] and 
connections with random matrix theory [38,44,181]), which could fit perfectly 
within the mainstream of this dissertation and will indeed have some bearing 
on the content of Chapter 
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The core of this section is the brief survey of spin CS models and Dunkl 
operators developed in Sections 11.61 and 11.71 While we introduce a wealth 
of relevant concepts, we do not include a single rigorous proof or statement: 
our presentation of these topics should be regarded as a large dry bone that 
will be fleshed in Chapters |2] and [3l 

We occasionally present some original results as we go along; these are 
mainly taken from Refs. [15,69,70]. 



1.2 Scalar CS models 

Let R C be a (crystallographic, possibly nonreduced) root system, which 
we assume to be irreducible without loss of generahty, and let i?+ C i? a 
subsystem of positive roots [118]. Given a. E R, we represent by 

Q! * X 

cra(x) = X - 2- — —a 

the reflection across the hyperplane orthogonal to a. The Weyl group cor- 
responding to R is denoted by 22J(-R) = (cTq, : ck G R+). 

Consider a function v G C°°(]R\{0}) of one of the following forms: 

(rational) 

v{x) = ^ sinh~^a; (hyperbolic) (1.1) 
sin~^ X (trigonometric). 

The (scalar) A^-body CS model of type R is defined by the Hamiltonian 

H,, = -A+ 9lv{a-^), (1.2) 

where 

N 



i=l 



and the coupling constants g'^ are assumed to take the same value for roots of 
the same length. In the rational case we shall frequently include a confining 
harmonic term, i.e.. 



^-^ (a ■ xr 
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with r = ||x||. If we define a unitary action K of W{R) on L^(M^) by 

K{(ra) = Ka ■■ f ^ f O (Ta, (1.3) 

it is clear that 2n(-R) is a symmetry group of the Hamiltonian (11. 2p . 

Remark 1.1. More generally [158], the function v in Eq. (11.11) can be taken to 
be a Weierstrass p function with suitable periods. Despite the relevance of 
this case and the abundant related literature (cf. e.g. [135,136] and references 
therein), we have preferred to omit elliptic potentials in our presentation of 
CS models. On the one hand, they do not appear anywhere in the sequel; 
on the other, their solvability and integrability properties are rather different 
from those of their limiting cases listed in Eq. (11.11) (with u set to zero). The 
conditions on the root system can be also relaxed somewhat [29,158]. 

The double pole that v possesses at imposes that any function ip in the 
domain of i^sc must vanish at the hyperplanes x - a = with a G which 
means that the particles cannot overtake one another. Hence one can assume 
that the configuration space for Hsc is the Weyl chamber 

C= {xGM^:a-x>0, VaGi?+} (1.4) 

in the rational and hyperbolic cases, and the Weyl alcove 

C = {x G : a ■ X > 0, Va G a^ax • x < vr} (1.5) 

in the trigonometric one. Here we denote by ctmax the highest root of the 
system R. 

Remark 1.2. The symmetry of under W{R) guarantees the equivalence of 
the eigenvalue problems for Hsc in and in C. Namely, for any eigenfunc- 
tion (f G L^(M^) of Hsc there exists an eigenfunction ip G L'^{C) of Hsc with 
the same eigenvalue and a one-dimensional representation x '■ 2n(i?) ^ {±1} 
of the Weyl group such that for any w G 2n(-R) satisfying wx G C one has 

(p{x.) = xiw) ^(wx) . 

Clearly, the inclusion map defines a one-to-one correspondence between the 
eigenfunctions of Hsc on LF'[C) and its eigenfunctions on L^(R^) with a given 
symmetry. 

We shall be mainly interested in the root systems A^v-i and BC^, which 
are respectively given by 



R= {ei-ej -.l^ij^j <:N} 
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and 

R={±ei, ±2ei, ±e, ± e,- : 1 ^ z < j ^ A^} , 

{ej : 1 ^ i ^ A^} being the canonical basis of M^. An appropriate choice 
of the couphng constants of the CS models associated with these systems 
suffices to exhaust all the possibilities arising from the non-exceptional Lie 
algebras. Physically, these systems are important because they are the only 
ones which allow for an arbitrary number A^ of particles, given by the rank of 
the algebra. For the sake of completeness and future reference, let us mention 
that set of positive roots i?+ and the highest root CKmax for the Ajy^i and 
BCn systems can be taken as 



R, 



{ei - ej : 1 ^ i < j ^ N} 



ei - Bat 



and 



R+ = {ei ± e^, ei, 2ei : 1 ^ i < j ^ N} 



Ctr 



2ei 



The An-1 and BCn scalar CS models (11. 2p with rational and trigonometric 
potentials can therefore be written as 



-f^sc = — A + Vsc , 

where the potential is in each case given by 



:i.6) 



a(a — 1) 2 2 



t * CI 



a(a — 11 



sin (xj — Xj) 



Vsc = a{a - 1) 
Vsc = a(a - 1) 



1 



1 1 



+ 



E 



(1.7a) 
(1.7b) 

!*^+-V, (1.7c) 



sin^(xj — Xj) sin^(xi + Xj] 



^ — ^ sm Xi ^ — ^ 



b'{b' - 1) 



COS^ Xi 



;i.7d) 



Here and in what follows the indices in all sums and products range from 1 
to A^ unless otherwise stated. 
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1.3 Integrability and solvability 

CS models are the paradigm of quantum integrable A^-body problems. The 
first satisfactory explanation of their remarkable algebraic properties was 
given by Olshanetsky and Perelomov [158], who constructed N commuting 
first integrals for the Hamiltonian (11.21) associated with any non-exceptional 
root system. A new proof of the integrability of these systems (which applies 
even to the case of non-crystallographic root systems) has been recently de- 
veloped by Bordner, Manton and Sasaki [31] using universal quantum Lax 
pairs, and previously by Heckman and Opdam [106, 107] (in the context of 
degenerate Hecke algebras) and by Oshima and Sekiguchi [162]. For sim- 
plicity, we shall henceforth assume that is a non-exceptional root system, 
typically A^^i or BCn- 

For these root systems, the models (11. 2p are exactly solvable in the 
sense that all their eigenvalues and eigenfunctions can be computed in closed 
form. (The solvability properties of CS models associated with exceptional 
Lie algebras are still uncertain, even though there have been some recent 
developments in this area [32].) The eigenfunctions are labeled by integer 
multiindices n = (ni, . . . , n^r) and factorize as 

V'„(x) =/i(x)P„(z), (1.8) 

where = ^/'o is the ground state function, z = (2:1, ... , zn) with 

Xi (rational) 
g2xj (hyperbolic) 

(trigonometric) , 



and Pn is a polynomial of degree 



|n| = E 



1.1 . (1.9) 



Physically, the multiindex n can be thought of as the quasimomentum of the 
particles [168,171], and the energy of the corresponding eigenfunction can be 
expressed neatly in terms of this quantity and the parameters appearing in 
the Hamiltonian. 

Example 1.3. The original Calogero model ( ll.7al) is obtained by setting 
R = An -I and choosing the rational potential with a harmonic term. Its 
eigenfunctions are labeled by a multiindex n e Nq and can be written as 
ipn = /^-Pn, where 



Mx) = e-t^^ n 



I o, 



i<j 
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is the ground state function and Pn is a homogeneous polynomial of de- 
gree |n|). These polynomials are known as generalized Hermite polynomi- 
als [137,200], since when = 1 they reduce to Hermite polynomials. Ex- 
plicit formulas can be consulted, e.g., in Refs. [34,103,125]. The energy of 
the eigenfunction ipn is given by 

En = 2u!\n\ + Eq , 

with Eo = ujN{a{N - 1) + 1). 

Example 1.4. If we keep the Atv-i root system and choose the trigonometric 
potential, we obtain the original Sutherland model. Its eigenfunctions can 
be written as ipni'^) = /^(x)Pn(z), where 

/i(x) = JJ|sin(xi-Xj)|", 

n G Z^, and Pn is a polynomial of degree |n| in z = (e^^^^, . . . , e^^^^). These 
polynomials are called Jack polynomials [123], and play a fundamental role 
in Heckman and Opdam's theory of hypergeometric functions for root sys- 
tems [104,109,159,160]. 

A hyperbolic version of the trigonometric Sutherland system can be read- 
ily obtained. Its spectrum, however, is not discrete, and in fact this model 
only has a finite number of eigenvalues. 

1.4 Orthogonal polynomials 

Whereas the theory of orthogonal polynomials in one variable is a classical 
branch of Mathematics, significant results on multivariate polynomials were 
not obtained until the second half of the XX century [61]. Since then, orthog- 
onal polynomials of several variables have become a major field of research in 
modern analysis and combinatorics [19,43,146,185,199]. As can be seen from 
Eq. fll.Sp . multivariate orthogonal polynomials appear as eigenfunctions of 
CS models, and the role that these models play in the latter theory is anal- 
ogous to that of the Laguerre, Hermite or Jacobi equations in the case of 
classical orthogonal polynomials in one variable. 

As a matter of fact, the analogy between the eigenfunctions of CS models 
and classical orthogonal polynomials is even closer. We showed in Exam- 
ple 11.31 that the eigenfunctions of the rational Aj^j^i CS model (the Hi- Jack 
polynomials) reduce to the Hermite polynomials in the case A^ = 1 and indeed 
coincide with their multivariate generalization [137]. It is also well known [14] 
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that the multivariate generahzations of the Laguerre [139] and Jacobi [138] 
polynomials are similarly recovered from the BCn rational and trigonometric 
CS models, respectively. Generally, CS models give rise to generalized hy- 
pergeometric functions associated with root systems, as extensively studied 
by Heckman and Opdam [104,109,159,160]. 

A recent development is some kind of converse statement to the latter 
comment. Hallnas and Langmann have constructed an algorithm that al- 
lows one to go from a family of orthogonal polynomials in one variable to 
an A^-body Hamiltonian whose eigenfunctions are a symmetric generaliza- 
tion of this family to variables. The A^-body models associated with the 
classical families of orthogonal polynomials are in fact previously known CS 
Hamiltonians. 

The fruitful connection between CS models and the theory of orthogonal 
polynomials has lead to a plethora of closed formulas for the eigenfunctions 
of these Hamiltonians (cf. [61,103,146] and references therein) . Furthermore, 
an exciting new ingredient of this theory are supersymmetric (spin) polyno- 
mials [52,53], related directly to the spin CS models which will be the main 
subject of the remaining chapters of this memoir. 

1.5 A classical detour 

It can be argued that the conceptually simplest approach to the exceptional 
solvability properties of CS models is through their classical counterparts, 
the Calogero-Moser (CM) models. Hence, we deem it important to present 
a brief account of the theory of CM models, which, apart from its academic 
interest, is essential in order to understand some recent developments on 
classical Haldane-Shastry spin chains [122]. 

CM models [157] are the Hamiltonian systems on whose Hamil- 
tonian function is obtained from Eq. (11.21) by substitution of the quantum 
kinetic energy (the Laplacian —A) for the classical one (the squared mo- 
mentum IIpP). The mathematical intricacies that lurk behind quantum 
integrability [45,69,134,201] are absent in the classical setting, the Liouville 
integrability of CM models being related to the existence of an invariant foli- 
ation by Lagrangian cylinders in phase space [2]. The Liouville integrability 
of CM models for any (generalized) root system stems from Bordner, Corri- 
gan and Sasaki's construction [30] of a universal Lax representation for these 
systems. Nonetheless, their integrability was also explicitly obtained [160] as 
a corollary in a series of papers by Heckman and Opdam ten years earlier 
by means of a detailed study of the corresponding CS models and a suitable 
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theorem on the classical limit of their quantum first integrals. An advantage 
of the first authors' approach is that the equations of motion of the sys- 
tem can be explicitly integrated for any root system (resp. for root systems 
associated with classical Lie algebras) when the potential is rational (resp. 
trigonometric or hyperbolic.) 

Most remarkably, the CM Hamiltonians with rational and hyper- 
bolic potentials are also maximally superintegrable for any value of N, 
i.e., they admit the maximum number 2N — 1 of functionally independent 
global first integrals. This property is quite exceptional: apart from these 
two models, the only known (natural) Hamiltonian systems which are max- 
imally superintegrable in arbitrary dimension are the geodesic fiow on the 
simply connected space forms, the (generalized) Kepler problem [73], the 
Smorodinsky-Winternitz system [74,90], their generalizations to the sim- 
ply connected spaces of constant sectional curvature (cf. [16] and references 
therein), the nonisotropic oscillator with rational frequencies, the nonperi- 
odic An Toda lattice [5,194], and a very recent example living on a manifold 
of nonconstant curvature [15], namely 

i7(x,p) = 



1 + r2 

This system can be understood as the intrinsic Smorodinsky-Winternitz 
model on an A^- dimensional generalization of the Darboux space of type 
III [131]. 

The exceptional properties of the CM models can be understood by means 
of (singular) Hamiltonian reduction [161]. This technique, independently de- 
veloped by Olshanetsky and Perelomov [155,156] and by Kazhdan, Konstant 
and Sternberg [130], was recently clarified and placed in a broader context by 
Feher and Pusztai [76-78]. The latter works provide some kind of symplectic 
analog of Harish-Chandra's theory of c-functions [72, 112]. 

We sketch the gist of the method for the simple case of the A^-i trigono- 
metric Sutherland modefl (cf. Example II. 4p . as developed by Kazhdan, Kon- 
stant and Sternberg [130]. Let us consider the geodesic motion on K = 
SU(A^), which is given by the free Hamiltonian 

Hoik,\) = \\\f 

on T*K, where we trivialize T*K = K x t* (with t = Lie{K)) by right 
translations and identify i* = t hj means of the inner product. Here the 



am indebted to Prof. L. Feher for valuable discussions on the state of the art of the 
projection method in CM models. 
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invariant symplectic form on T*K is = d(A, {dk)k The Hamiltonian 
action 

{k,\) ^ {hkh-\hXh-^) , heK (1.10) 
is a symmetry of Hq and gives rise to a momentum map ijj : T*K — > i defined 

by 

Now let us look at the (singular) Marsden-Weinstein reduced phase spaces 

where denotes the isotropy group of /i G 6. (Generally speaking, these 
symplectic spaces are not manifolds, but stratified sets whose strata are dif- 
ferentiable manifolds.) The easiest nontrivial choice for is 

/i = i(W-iv-V), 

u G being an arbitrary vector of norm g. In this case, the reduced 
phase space actually turns out to be = T*C = C x Cartan(i^') = C x 
R-'^, C being a Weyl alcove, and the reduced Hamiltonian is given by the 
trigonometric Sutherland model, with (x, p) G C x M^. The integrability 
and solvability properties of the latter system now follow from those of the 
geodesic flow on K. As a matter of fact [77] , these ideas can be also applied 
to obtain a Lax representation. These ideas apply to any compact, connected 
simple Lie group K, but actually it can be proved that a scalar CM model 
is only obtained in the case that we have reviewed, whereas for a general K 
a spin CM model arises. 

The above construction can be reinterpreted as follows. Let G = K x K 
and consider its involution Q{ki,k2) = (^2,^1)- If = {{k,k) : k G K} 
and G^ (resp. G^) denotes its right (resp. left) action on G, it is clear that 
geodesic motion on G enjoys a G^ x G^ symmetry. If we reduce T*G using 
the G^ symmetry at the zero value of the momentum map, we arrive precisely 
at the symplectic system {T*K, u, H) analyzed in the preceding paragraph. 
Moreover, the G^ symmetry still survives, and its action on ~ G/G^ 
in fact coincides with the conjugation action fll.lOp of K on itself. In this 
framework, hyperbolic CS models (possibly with spin) are obtained when 
considering the complexification G oi K given by g = 6 + it, B being a 
Cartan involution of G. The hyperbolic Sutherland model is obtained when 
K = SU(A^) ~ G+, so that G = SL(iV, C), G/G^ ~ exp(it) and = T*G 
for the above prescription for /i, C now being a Weyl chamber. 

The above results also hold true, mutatis mutandis, if we take G to be any 
noncompact simple Lie group with finite center, G a Cartan involution, and 
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= {g E G : Q{g) = g}, so that it is a maximal compact subgroup. Any 
element g E G can be uniquely decomposed as g = g+g-, where g± G G±. and 

= {g E G : Q{g) = g^^} is diffeomorphic to 0_ via the exponential map. 
In this framework, one usually obtains spin CM models on T*G, and further 
generalizations have been also considered. Moreover, the quantum analogs of 
these results are being actively pursued by the same authors, and connections 
with flat symmetric spaces have also been studied recently [6,115]. 

The related topic of the reduction of CS models by discrete symmetries 
has also been considered in the Physics literature [167,170], both in the spin 
and spinless cases. A wide panoply of CS models can be obtained from the 
original Calogero model through this procedure, partially explaining why 
this system (with or without the harmonic confining potential) is sometimes 
regarded as the mother of all CM systems [171]. 

Interestingly, some purely quantum mechanical techniques in CS models 
were originally suggested by a careful analysis of their classical analogs. For 
instance, Polychronakos's rediscovery of Dunkl operators [164] was motivated 
by a careful derivation of the scattering Calogero model by reduction of free 
motion. Moreover, a short calculation allows one to prove that the only 
asymptotic effect of the interaction with the potential fll.7al) is an effective 
reshuffling of the particles. This property carries over to the quantum case, 
where it acquires an even greater relevance through the concept of identical 
particles. It is therefore amusing that some "physical" effects of the CS mod- 
els were first observed at a classical level but only become really fundamental 
in the quantum realm. For instance, the ideas that we have just sketched lie 
at the heart of the applications of CS models to generalized statistics in one 
dimension. 

Another fruitful connection between CM and CS models, which actually 
holds for a wide variety of one-dimensional quantum systems, appears when 
considering the exact bosonization of one such model. This technique 
essentially consists [27,195,203] in describing a fermionic system by means 
of a bosonic field theory. More precisely, we require that the algebra of 
bosonic operators be irreducible on the fermionic Hilbert space and that the 
bosonic field theory reproduce the exact eigenvalues of the system up to non- 
perturbative effects. It has been recently discovered [70] that an exact one- 
dimensional bosonization can be constructed via Weyl quantization of the 
classical picture, thus obtaining a new realization of the W^o algebra. This 
result is expected to shed new light on higher dimensional bosonization, but 
the difficulties encountered in this case are considerable. 
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1.6 Dunkl operators 

The first attempt to provide a conceptual explanation of tlie exact solvability 
of CS models was made by Olshanetsky and Perelomov [158] . In the rational 
(with ci; = 0) and trigonometric cases, they showed that there exist some 
particular values of the coupling constants for which the Hamiltonian fll.2p is 
equivalent to the radial part of the Laplacian on a symmetric space associated 
with the root system R. The radial parts of the Laplacians on these spaces 
have been thoroughly studied [112], and quite explicit formulas are available. 
The general rational {u ^ 0) and hyperbolic cases could also be dealt with 
using appropriate modifications. The main drawback of Olshanetsky and 
Perelomov's approach was that for the BC^ case and arbitrary coupling 
constants one did not obtain the radial Laplacian. This, however, has been 
recently amended by Feher and Pusztai (although so far only the classical 
case has been carried out in detail [78]). 

Concerning the connection between CS models and the theory of orthog- 
onal polynomials, a major breakthrough was the introduction of Dunkl (or 
exchange) operators [59] 

jr' = d.^-y] (1 - K^) , (1.11) 

^ — ^ Q. ■ X 

to study spherical harmonics associated with measures invariant under the 
Weyl group W{R). Dunkl operators have become an essential tool in har- 
monic analysis, combinatorics and CS models [107, 175]. The differential- 
difference operators (11.111) have two key properties: they are homogeneous of 
degree —1, i.e., their action on a polynomial lowers its degree by one, and are 
mutually commutative. In this sense, one can think of the Dunkl operators 
as deformations of the usual partial derivatives that preserve most of their 
analjd;ic properties. In fact, one can construct a Dunkl transform [46, 47] 
using the spectral resolution of the anti-self-adjoint operator J'''^*, and the 
resulting properties are fairly similar to those of the usual Fourier transform. 

In the Physics community, Dunkl operators were rediscovered by Poly- 
chronakos [164] within the context of CS models. To explain the essentials 
of the connection between these objects, let us restrict ourselves to the case 
v{x) = x^"^ and define 

/i(x) = n I" • 

A simple link can now be established through the equation 

fiH,,fi^' = -Y,{Jn' + Eo, (1.12) 
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where Eq is a constant and the constants ga are given by 

This equation is known to hold true [92] on the space L^(R^)^^^) of 2U(-R)- 
invariant square- integrable functions on M^, which is canonically isometric 
to our (reduced) Hilbert space L^(C). Hence Dunkl operators appear in the 
rational {uj = 0) case (the scattering Calogero model) as first integrals of 
the system, and have the physical interpretation of (conserved) asymptotic 
momenta. Furthermore, Eq. (11.121) partially accounts for the relevance of CS 
models in harmonic analysis. External harmonic potentials can be introduced 
in this picture using either deformed destruction operators Oj = —i{Jl^^ + 
uxi) [164] or auxiliary differential operators [81]. 

In the trigonometric case one can define an analogous set of commuting 
differential- difference operators by 

The relation between the Sutherland system of type R and the sum of squares 
of these Dunkl operators goes along the lines of the rational case (11.121) . with 
gauge factor 

^(x) = Y[ |sin(a-x)P" . (1.14) 

It can be readily verified that the operators (11.131) leave invariant the subspace 
of polynomials in z = (e^^^^, . . . ,e^'^^) of degree not higher than k, for all 
k G No- 

As a matter of fact, the trigonometric Dunkl operators (I1.13P were in- 
troduced by Cherednik [41,42], while a closely related set of differential- 
difference operators has been previously considered by Heckman [105]. Sev- 
eral analogous families of operators are currently known, both in the rational 
and trigonometric cases. For simplicity, we shall understand by (generalized) 
Dunkl operators of type R [81, 107, 175] any set of first order differential- 
difference operators Jj (z = 1, . . . , A^) in variables z satisfying: 

(i) The operators Jj map the space of polynomials C[z] into itself and leave 
invariant finite-dimensional subspaces A^" C C[z]. 

(ii) The operators { Jj, Ka : 1 ^ i ^ N, a E R+} span a degenerate Hecke 
algebra. 
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Degenerate (or graded) Hecke algebras were defined independently by Drin- 
feld [58] and Lusztig [144]. For the sake of completeness, let us recall 
that the degenerate Hecke algebra (H, *) associated with the system of pos- 
itive roots -R+ C V and the parameters {aa)ai=R+ is the linear space H = 
S(y'^) ® C[2IJ(-R)] endowed with the unique graded algebra structure satis- 
fying 

(i) S{V^)®1 and mC[W{R)] are subalgebras, so that S{V^) and C[W{R)] 
can be identified with their images in H. 

(ii) If CK is a simple root, then <7a*"v — {(Ta'v)*(Ja = — OqCk-v for all v G V^. 

Here R+ is assumed to span V, Siy^) denotes the algebra of symmetric 
tensors in the complex linear space V^, and the grading in EI is in fact 
inherited from that of S{V'^). Clearly one can identify T *w = T ^w, where 
T e SiV^) and w E C[W{R)]. 

As in Eq. (11.121) . Dunkl operators happen to give rise to (generalized) 
CS models through the addition of auxiliary operators, a similarity trans- 
formation and a global change of variables. Full details can be consulted in 
Refs. [81,82] and in Chapter El 

1.7 Spin CS models 

We shall comment now on another benefit granted by the Dunkl operator 
formalism that is essential for all the forthcoming sections: the definition of 
spin CS models. Since the study of (generalizations of) such models could 
be rightly considered to be the core of this Ph. D. dissertation, we shall have 
many opportunities in the following chapters to elaborate on this topic. Our 
goal in this section is to define the basic operators and present the funda- 
mentals of the theory without providing any details on the computational 
difficulties involved. This and other aspects shall be covered in detail in 
Chapter El 

For the sake of simplicity, we shall restrict ourselves to the non-exceptional 
root systems An and BCn, which are the only ones that have been considered 
in the literature up to now. (In fact, we shall make use of the root system 
An-1 embedded in M.^ , not of A^. Nevertheless, we shall hereafter drop the 
subscript —1 and regard the subscript as an abstract index.) 

Let S ^ 0^^+^)"" be the Hilbert space of the internal degrees of freedom 
of particles of spin M e iN, and let 

B^ = {\s) = \s,,...,SN):s,,e {-M, -M + 1, . . . , M}} (1.15) 
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be the canonical basis of S. Let us define the spin permutation and reversal 
operators by 

Sij\si, . . . ,Si,. . . ,Sj . . .sn) = \si, ...,Sj,...,Si...SN), (1.16a) 
Si\si, . . . ,Si, . . . sn) = \si, -Si, ... Sat) , (1.16b) 
Sij\s) = SiSjSij\s) . (1.16c) 

One can easily verify that they satisfy the algebraic relations 



i~>ij — -L; iJijiJjk — >->iki~>ij — i~>jki~>ikj >->ij>->kl — >->kl>J', 



q'2 -| CO Q C C Q CO O O O O 

where the indices k, I take distinct values in the range 1, . . . ,N. 

The one-dimensional representations x '■ 2n(-R) — {±1} of the Weyl 
group W{R) are given by [118] 



e = ±1 



in the An- case and by 



0"p 



e, ■ ' - , e, e' = ±1 



in the BCn case. For each of these one-dimensional representations x? oiie 
can define the spin CS model with scalar counterpart fll.6p as 



-A + K 



where the spin potentials Vy are given by 



, sin^fxj — Xj) ' 



E 
E 



a(a — eSii) ^ a(a — eSj, 



(x^ "^j)^ i^Xi ~\- Xj^'^ 



a(a — eS*,-. 



a(a — eSi. 



sin^(xj — Xj) sin^(xi + Xj) 



E 



b{b - e'Sj 
sin^ Xi 



+ E 



b'{b' - e'Si) 



COS^ Xj 



(1.17) 
(1.18a) 

;i.i8b) 

('1.18c) 



;i.i8d) 
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We shall frequently omit the subscript if there is no risk of confusion. The 
extension to the hyperbohc case is straightforward. 

Let us introduce the coordinate permutation and sign reversal operators 

{Kijf){xi, . . . ,Xi, . . . ,Xj, . . .,xn) = f{xi, . . . ,Xj, . . . ,Xi, . . .,xn) , (1.19a) 
{Kif){xi, ...,Xi,.. .,xn) = f{xi, -Xi, . . .,xn) , (1.19b) 
ki,f = K,KjK,jf, (1.19c) 

which are the reflection operators associated with the root systems Ajy and 
BCn- One obviously has the group isomorphisms (Kij) ^ (Sij) ~ W{Aiy_i) 
and {Ki,Kij) ^ {Si, Sij) ^ W{BCn). We shall denote by i^, e and W re- 
spectively the realizations of 2n(-R) generated by the operators (11.191) . (11.161) . 
and 

Uij=Ki,Sij, Ui = KiSi. (1.20) 

When notationally convenient, we shall regard these realizations as unitary 
representations 

K : w ^ e End{L^) , (1.21a) 
S -.w^ Sy,e End(S) , (1.21b) 
n : w ^ G End(L=^ ® E) , (1.21c) 

with w G W{R). Moreover, we shall not distinguish between a character of 
W{R) and its isomorphic extension to 2IJ, 6 or 

The key property of the spin models (11.181) is their connection with Dunkl 
operators. This connection is implemented via the star mapping 
il((5) associated with the one- dimensional representation x, which is the 
antihomomorphism defined by 

Kij ^ (Ki,):^, = tSij , (1.22a) 
Ki ^ {Ki)l, = e'Si . (1.22b) 

Here il(g) stands for the (real) universal enveloping algebra of g, and the 
second equation fll.22bp is to be omitted in the A^ case. We shall also 
consider the extension of (11.221) to D^^^D^ & given by 

(D ® K):, = D® (K):, . 

Here D is the space of (scalar) smooth differential operators, and D and 
K belong to D and K respectively. It is easy to show that the spin CS 
models (I1.17P are obtained by applying the star mapping associated with x 
to the exchange Hamiltonian 

H=-A+V, (1.23) 
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with 




(1.24b) 



(1.24a) 




a(a — Kij) a{a — Kij) 



a{a — Kij) a(a — Kij) 




h{h - Kj) 



(1.24c) 





+5: 




(1.24d) 



Since these are essentially the operators obtained from the sum of the squares 
of the Dunkl operators (11.111) and (11.131) through a gauge transformation and 
a change of variables, the eigenvalue problem for the spin Hamiltonian (I1.17P 
can be effectively reduced to the analogous problem for a quadratic polyno- 
mial in the Dunkl operators plus, perhaps, some simple auxiliary operators. 

The connection between CS models and Dunkl operators can also be used 
to define generalizations of the above solvable Hamiltonians (with or without 
spin). Such an operator H is called of CS type. Its spectrum cannot be 
generally computed by algebraic methods, but by construction [81,82] H 
remains partially solvable in the sense that it leaves invariant some known 
finite-dimensional subspace of L^(M^) ® S; in this case, H is said to be 
quasi- exactly solvable (QES), and the eigenvalues that one can obtain 
from the latter subspaces are called algebraic. 

Remark 1.5. A particularly interesting situation [196,197] arises when one 
can find explicitly an infinite flag (i.e., a sequence A^i C S ' ' ' C 
L^(M^) ® S, where Ain is finite-dimensional) invariant under H. In this 
case, the corresponding Hamiltonian H is sometimes said to be exactly solv- 
able in the sense of Turhiner [196]. This expression is certainly misleading, 
the latter condition does not imply by any means that algebraic eigenvalues 
cover the whole point spectrum of H . In this dissertation we shall be inter- 
ested in this stronger type of quasi-exact solvability, and in fact all the QES 
Hamiltonians that we shall encounter in Chapter [3] will indeed preserve an 
infinite-dimensional invariant flag. 
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1.8 Solvable spin chains 

Solvable spin chains have enjoyed a growing popularity in the last few years, 
due in part to their novel applications to SUSY Yang-Mills and string the- 
ories [20,21,89,94,150,173]. The role of CS spin models in the study of 
solvable spin chains was unveiled by Polychronakos [164, 166] in a couple of 
insightful papers from the mid 90s which elucidated the connection between 
the original Sutherland model and the recently discovered Haldane-Shastry 
(HS) chain [100,178]. We shall conclude this chapter with an overview of 
the essentials of the CS/HS connection and the definition of some celebrated 
spin chains that shall be frequently referred to in forthcoming chapters. 

The Heisenberg chain [110], which was born as an attempt to model 
ferromagnetic materials, describes particles on a lattice with isotropic near- 
neighbor interactions independent of the site. For particles of spin M G |N, 
the Hamiltonian of the system is customarily written as 

HHe = $^S,-S,+i, (1.25) 

i 

where Sj is the (vector) spin operator of the particle i, Sn+i = Si and 

4M(M+1) 
j ■ Dj — Oj Dj . 

a=l 

For each fixed particle label i, we denote by (a = 1, . . . ,4:M{M + 1)) a 
basis of the fundamental representation of su(2M + 1), which we assume to 
be orthonormal with respect to the bilinear form {S,S') = 2tr(S'S"). It is 
straightforward to write the Hamiltonian (I1.25P in terms of the spin exchange 
operators fll.l6ap using that 

S^, = 2 S, ■ S,- + — . (1.26) 

^ * J 2M+1 ^ ' 

Note that for spin 1/2 particles we have Si = ^{a},af,af), where aj is the 
j-th Pauli matrix acting on the spin space of the i-th particle. 

It is well known [25,50, 117] that the spin 1/2 Heisenberg chain can be 
solved exactly using the Bethe ansatz. Several (partially) solvable gener- 
alizations of the Heisenberg chain (11.251) with short-range interactions (at 
most between next-to-nearest neighbors) have been subsequently proposed 
in the literature. These include, in particular, the family of chains with 
arbitrary spin and nearest-neighbors interactions polynomial in Sj ■ Sj+i of 
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Refs. [12,192], as well as several models whose ground state can be written 
in terms of "valence bonds" [3, 147, 148]. 

The Haldane-Shastry model [100, 178] is the prime example of solv- 
able spin chain with long-range interactions. It describes an arrangement of 
N equally spaced spins on a circle which interact pairwise with a strength 
proportional to the inverse square of their chord distance. The Hamiltonian 
of the HS chain can be written as 



where the positions of the sites are given by 



This system quickly attracted considerable attention [86,101,120,129,188] 
due to its remarkable algebraic properties and its connection with conformal 
field theory [55] and the symplectic Dyson model [62]. The original interest 
of this model, however, lies in the fact that Gutzwiller's variational wave 
function for the Hubbard model [91,97,98], which also coincides with the 
one-dimensional version of the resonating valence bond state introduced by 
Anderson [7] and with the one-dimensional version of the Kalmeyer-Laughlin 
state [126], becomes an exact eigenfunction of the HS chain as the strength 
of the on-site interaction tends to infinity. 

In the early nineties the integrability of the Hamiltonian fll.271) was rig- 
orously proved and its Yangian symmetry [87, 102] was established. Not 
surprisingly, the Yangian algebra [57] also appears as the symmetry alge- 
bra [9,23,113] of the Sutherland model ( 11.7b[) . Although the obvious re- 
lationship between this model and the chain fll.27p was already remarked 
by Shastry in his original paper, their connection was first made explicit by 
Polychronakos through the so-called freezing trick [165]. This technique 
formally allows one to recover the HS Hamiltonian from the Sutherland spin 
model (11.18bl) (with e = 1) in the large coupling constant limit, which is 
clearly tantamount to the classical limit /i | 0. A rigorous proof of the 
classical limit and some extensions and improvements will be presented in 
Chapters [2] and H 

A further application of the freezing trick enabled Polychronakos to obtain 
the integrable spin chain with long-range interactions given by the Hamilto- 
nian 




(1.27) 




(1.28) 
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where the sites < ■ ■ ■ < ^at satisfy the equation 



6 = E 



1 



(1.29) 



- 



This system is usually referred to as the Polychronakos-Frahm (PF) 
chain, and relates to the spin Calogero model fll.lSal) as the HS chain does 
to the Sutherland model fll.lSbp . 

As remarked by Frahm [88] , the sites of the PF chain are no longer equally 
spaced, but Eq. f ll.29p shows that they are given by the zeros of the A^-th 
Hermite polynomial [189]. In particular [36], this implies that the density of 
sites for large is given by 



Remark 1.6. The CS/HS correspondence casts the original model of Haldane 
and Shastry into a much wider framework. Consequently, we shall say that 
a spin chain is of HS type if it can be derived from a spin Schrodinger 
operator of CS type via the freezing limit. In particular, both the HS and 
the PF chains are of HS type. 

Interestingly, we shall see in Chapters [2] and H] (cf. also [84,166]) that 
the chains of HS type, featuring long-range position-dependent interactions, 
seem to enjoy stronger solvability properties than those of Heisenberg type, 
characterized by the short range and position independence of the interac- 
tions. 



pP^(a;) = 4-y2iV^. 



(1.30) 
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Chapter 2 

HS spin chains of BCjsf type 



2.1 Introduction 

In the previous chapter we saw that the connection between An CS models 
and spin chains of HS type has been thoroughly studied. In contrast, BCn 
spin models [39,60,82,122,205] had received comparatively little attention. 
Two main reasons were responsible for this fact. On one hand, BCn HS 
chains depend nontrivially on free parameters (one in the rational case and 
two in the trigonometric one). On the other hand, the sites of the BCn 
chain are not equally spaced and no explicit description of their positions is 
available. 

Let us give a very brief review of the existing results on these models. The 
integrability of the rational BCn HS chain was established by Yamamoto 
and Tsuchiya [206] using Dunkl operators, but its spectrum has not been 
computed so far. The HS trigonometric spin chain of BCn type was discussed 
by Bernard, Pasquier and Serban [24], but only for spin 1/2 and with the 
assumption that the sites are equally spaced, which restricts the space of 
free parameters in the model to just three particular values. Finkel et al. [83] 
discussed the integrability of the hyperbolic HS chain of BCn type, but did 
not examine its spectrum either. 

The main result of this chapter is the exact computation of the spec- 
trum of the BCn HS chain. In fact, four chains associated with the BCn 
Sutherland model are considered, two of which are essentially different. This 
computation is based on a rigorous study of Polychronakos's "freezing trick", 
which is in fact promoted to a "freezing lemma" in Section 12.51 

The organization of this chapter is as follows. In Sections 12.21 to 12.41 we 
present a detailed study of the BCn Sutherland spin model and prove its 
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exact solvability using Dunkl operators. In Section 12.51 we introduce the 
BCn HS chains and compute their spectrum, whose statistical properties we 
analyze carefully in the following section. Finally, in Section 12.71 we consider 
the addition of a constant magnetic field to this picture. 

The material presented in this chapter is based on Refs. [63-65]. 

2.2 The spin BCn Sutherland model 

In this section we shall study in some detail the spin BCn CS models associ- 
ated with the one-dimensional representations Xee' of W{BCn)- Physically, 
the parameter e determines the ferromagnetic (e = 1) or antiferromagnetic 
(e = — 1) behavior of the associated spin chains. 

Let us recall that the potential Ke' of the spin BCjq Sutherland model 
was given in Eq. (jl.lSdp . For reasons that will be apparent by the end of this 
chapter, it is convenient to define the BCn Sutherland Hamiltonian by 
the action of the differential operator 

H,, = -A + Ke' (2.1) 
on the "reduced" Hilbert space 

n = L^{C)®^. (2.2) 
Here the Weyl alcove fll.5p reads 

C = {x : < a;i < • ■ ■ < Xiv < 7r/2} . 

For the sake of completeness, we shall prove rigorously that this choice of 
domain is perfectly admissible from a mathematical point of view. To this 
end, let us define the C^ function yU as in Eq. fll.l4p . so that 

/i(x) = JJ^ I sin(xi — Xj) sin(xj + x^)]" JJ^ | sinxi|^ JJ^ | cosXj|^ . (2.3) 

i<j i i 

Here we are assuming that a,b,b' > 1/2. It can be readily verified that 

where the scalar Hamiltonian Hsc is defined on L?'{C) by Eqs. fll.6l) - (]1.7dp 

and 

Eq = N[l{N - 1)(2A^ - l)a^ + 2{N - l)a{b + b') + (6 + b'f] . (2.4) 
Since /i does not vanish on C, fi is the ground state function of ifgc [143]. 
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Proposition 2.1. If a,b,b' > 2, H^^i is essentially self-adjoint onV^J], 
where 

Proof. A theorem of Komori and Takemura [132] shows that "D is a core for 
the operator Hsc ( 11.7d[) . Since H^^' — Hsc is i^sc-bounded with relative bound 

max{(a - iy\ {h - {h' - l)-^} , 

the Kato-Relhch theorem [172] shows that H^^i is essentially self-adjoint on 
E. □ 

Remark 2.2. Although its proof is simple (in fact, a direct proof could have 
been provided without much additional effort, and the conditions on a, 6, h' 
could have be easily relaxed), the inverse-square singularity prevents us from 
obtaining Proposition 12. II from the usual theorems on self-adjointness [184]. 
Interestingly, the proof of this result is essentially algebraic, and leans on our 
precise knowledge of the eigenfunctions of the BC]^ scalar Sutherland model. 

Although it has some applications to CS models [17,18,79], hereafter we 
shall not touch upon this topic, and by "self-adjoint" we shall always mean 
"formally self-adjoint". 

Some words on the choice of Hilbert space are in order. Consider the 
orthogonal projectors 

^"' = ^E^«'(n)n, (2.5a) 

A. = ^ 5^ x.(n)n, (2.5b) 

which we shall call symmetrizers , and let 

= {x : -7r/2 < Xi < n/2} 

be the fiat iV-torus. Observe that the angle coordinates of the torus are given 
by 2x, so that every function in C°(T^) will be vr-periodic in Xj. Define the 
"physical" Hilbert spaces of distinguishable particles as 

Ho = L2(T^) ®S. 

It is a standard result in representation theory [183] that Tio can be decom- 
posed into the direct sum 

Ho= n,e' (2.6) 
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of the auxiliary spaces 



which are given by the image of the symmetrizers fl2.5p . We shall also be 
interested in the subspaces 



i.e., the fermionic (e = —1) and bosonic (e = +1) sectors of Tio. 

It is technically convenient to define other two self-adjoint operators de- 
fined by the differential operator (12.11) . We shall denote by iJ*^*^ (resp. T^^') 
the Hamiltonian defined by the action of (12. ip on Hq (resp. Hee')- In the 
next proposition we summarize some elementary results relating the "physi- 
cal" Hamiltonians to the "reduced" ones. 

Proposition 2.3. The following statements hold: 
(i) 22J maps Hee' into itself. 

(a) W is a symmetry group of the physical Hamiltonian if"'. 

(Hi) The decompositions (12. 6p and i\2.7\} are invariant under the physical 
Hamiltonians. 

(iv) H^^'lt^^^, and H^^i are isospectral, and the corresponding eigenfunctions 
eH and $ G 7"^^^' are related by 



where w G W{BCiy) is such that G C . 

Proof. The first part of the lemma follows most easily from Lemma 12.51 in 
the next section. To prove the remaining points it is convenient to keep the 
abstract root formulation, with = \\cx.\\'^aa{aa — 1), and use the represen- 
tation (]1.21bp . Then one can write 



Tie = 'He.+ ® T^t- ) 



(2.7) 



$(x) = Xm'{w)Sw'^{wx.) , 



(2.8) 




(2.9) 



(2.10) 
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so that [i/" , n^] = 0. Since the decompositions (12.61) and (12. 7p are associated 
with a symmetry group of W^*^ , they are invariant under this operator. 

Each eigenfunction <I> of W^^' with energy E belongs to some subspace 7^^^/ 
by the symmetry. The potential diverging quadratically at dC C T^, $ 
must vanish at dC, and therefore \1/ = $|c must be an solution to the 
PDE 

(//«/-E)^ = inC, 
'^\sc = 0. 

Hence is an eigenfunction of the self-adjoint operator H^^i and $ can be 
recovered from \1/ through Eq. (12.81) . □ 

Corollary 2.4. spec(Tee') = spec(ifee')- 



2.3 BCn Dunkl operators 

In this section we shall make use of the Dunkl operator formalism to compute 
the spectrum of the spin model H^^i in closed form. Following Refs. [83, 
205] and Eq. (ll.lSp . let us consider the following set of commutative Dunkl 
operators 

Ji = id^. + a 1^(1 — i cot(xj — Xj)) Kij + (1 — i cot(xj + Xj)) Kij 

+ [b{l-icotXi) + b' {l + iteiiaXi)]Ki-2a^Kij. (2.11) 

j<i 

The following commutation relations with the exchange operators Ki,Kij 
can be readily checked: 

(2aiKik - K,k)Kij , iii<k<j 
[Kij, Jk\ = < (i.liaj 
[ , otherwise 

Ji - JjK.j = 2a(^l+J2 K^jKi^ , (2. 12b) 

i<i<j 

KijJ, - JiK.j = -2a (l + 5Z Ki,K,i^ , (2.12c) 

i<i<j 

[Ki, J,] = 2aKijiKi - Kj) , [K„ J,] = , (2.12d) 

{Ki, Ji} = 2(6 + b') + 2aJ2 Kii{Ki + Ki) , (2.12e) 

l>i 
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Here we assume that i < j and k are distinct indices in the range 1, . . . , A^. 

As an apphcation of the Dunkl operator formahsm, let us show the com- 
plete integrability of the spin Hamiltonian H^^i. The proof is based on the 
following lemma, which has been essentially taken from Ref. [83]. The def- 
initions of the spaces and (5 and of the star mapping were given in 
Section 11.71 

Lemma 2.5. Let A,B e D ® K and T e 2) (g) 6. Then the following 
statements hold: 

(i) AA,,, = 4VA«'. 

(ii) IfTA,,f = 0, then T = 0. 

(in) If B commutes withA^^i, then {AB)l^, = Al^,B*^,. 

(iv) IfA,B commute with A^^i , then [A,B]*^, = [A*^,,B*^,]. 

(v) If A commutes with K^j for some w G W{BCn), then A*^, commutes 
with 

Proof. (i) Let us write 

A = J2DaK^, (2.13) 

where a takes values in W^BCn). By definition, K^A^^i = Xee'i^r) SaA^^', 
so 

AA,,^ = J2xeei(y)D^S^A,,, = • 

(7 

(ii) Write T = E|n|<fc,a 5a,n(x)9"5^, where n G and 9" = d^l ■ ■ -Sjf^. 
The statement follows by choosing a finite set of functions {^Pn,a} C Het' 
such that the matrix (Sad^ip^'yi^)) is nonsingular. 

(iii) It follows from (i) and (ii) using that 

ABA,,, = AA,,,B = Al^,A,,,B = Al^,B*^,A,,, . 

(iv) It is a direct consequence of (iii). 

(v) Let us denote by Ad^ the automorphism of W{BCn) given by a i— *• 
waw. By hypothesis, and using the notation fl2.13p . 

A = K^AK^ = K^{D,)K a , 
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so Kyj{D(j) = Dj^d^a for each a. Hence 

<J 
a 

as we wanted to prove. 



□ 



Let us define the operators 

h = Y.^'" (2-14) 

i 

acting on L?'[C\ It can be readily verified that I\ coincides with the exchange 
Hamiltonian ( 11.24d[) . so that the spin Hamiltonian ( ll.lSdl) can be recovered 
as a differential operator from the relation H^^i = (h)*^/- Since the Dunkl 
operators fl2.11l) are self-adjoint and commute among them, so do the positive 
operators Ip. 

Proposition 2.6. The operators f l2.14p are ^-invariant. 

Proof. Let us prove that the elementary permutation -ft'j^i+i commutes with 
Ip. Eq. fl2.12al) shows that [i^j^j+i, Jj] = for j i,i + 1. Since 

2p-l 

r=0 
2p-l 

r=0 

it stems that [i^'j^j+i, ^ jf^] = 0. also commutes with Ip by virtue of the 
relations 

KNJi = JiKN, ifi < N , 
KnJ'n = —JnKnJn + 2(6 + h').JA! = JnKn ■ 

Here we have used that {K^, Jn} = 2(6 + b'). Now it immediately follows 
that 

[Ki, Ip] = KiN[KiNKi, Ip] = KiN[KNKiN, Ip] = 
for alH = 1,...,A^. □ 
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Lemma 2.7. The operator {Ip)*^i is self-adjoint. 

Proof. Proposition 12.61 shows that Ip commutes with A^e'. Using Lemma [2.51 
and the self-adjointness of Ip and A^^', one can see that 



Thus we have proved the following 

Theorem 2.8. The self-adjoint operators {Ip)*^i (p Efi) form a commuting 
family which includes the Hamiltonian H^^i = {Ii)*^>. 

2.4 Spectrum of the dynamical model 

There are two steps in the computation of the spectrum of H^^/. First, a key 
observation is that the gauge Hamiltonian is given by 



of finite-dimensional subspaces. The action of H is triangular in a suit- 
ably chosen basis and the eigenvalues can be read off the diagonal elements. 
Second, a careful analysis allows us to extend these results to the spin Hamil- 
tonian H^^/. 

Let us begin with some definitions. Let us consider the flag 



so [(/p)*^/]^ = {Ip)l^' as a consequence of Lemma [231 



□ 




(2.15) 



and that Jj leaves invariant an infinite flag 



C C ■ ■ ■ 



where 



n" = (A : |n| ^ k) . 



(2.16) 



Here n G Z 



n| is defined as in Eq. (11. 9p and 



/n(x) = /i(x)e' 



,2inx 



(2.17) 



with the gauge factor /i given by Eq. (12.31) and 



n ■ X = 
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Using Fourier analysis it is easy to see that the countable union IJfcLi is 
dense in L^(T'^), and that one can obtain a Hilbert basis E of this space by 
applying the Gram-Schmidt procedure to the functions fl2.17p . 

Given a multiindex n G Z^, we shall define the non-negative, nonincreas- 
ing multiindex 

[n] = {\n-„(i)l \n^{N)\) , 
where the permutation tt G is chosen so that 

The set of non-negative, nonincreasing multiindices will be denoted by [Z^]. 
We shall order by defining n -< m if there exists a number i ^ 1 such 
that [n]j = [m]j for all j < i and [n]j < [m]j, where [n\j denotes the jth 
component of [n]. This partial order can be naturally extended to {/n} and 
hence to B. We shall assume that the ordering of B is compatible with this 
partial order. 

In the following lemma we show that the order relations do in fact hold 
on orbits of ^. 

Lemma 2.9. For each w G W{BCiy), K^f^ = fwn and [wn] = [n]. In 
particular, there exists as element K & R such that Kfn = f[n] ■ 

Proof. The first part is trivial from the fact that w is a reflection, and thus 
n ■ wx. = wn ■ X. The second statement holds because W{BCn) is generated 
by Kij and Ki, i.e., by permutations and sign reversals. □ 

Now let us introduce some more notation. Given n G [Z^] and k G No, 
we define 

#(A;) = n) = card{i : = k}, 
e{k) = i{k, n) = min {{i : m = k} U {+oo}). 

Now one can characterize the action of the Dunkl operators (12. lip on the 
functions f l2.17p as follows. 

Proposition 2.10. For each n G [Z^], the action of Ji on the function /„ 
is given by 

Jifn ^i,nfn ~l" ^ ^ Cj j^/m , 

where 

^ _ (2n, + b + b' + 2a{N + i + l-i^{ni) -2i{ni)), if Ui > , 
\^b - b' + 2a{i - N) , if Ui = , 
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and G M. 

Proof. Let Zi = yij = ZiZj'^ y Wij = ZiZj. A tedious but straightforward 
computation shows that 



Jifn 2/n 



+ aiN-l) + Ub + b')+aJ2i-^^+ 



n. 

2 - 



Vij - 1 Wij - 1 



It can be readily verified that all the multiindices appearing in this formula 
are smaller than or equal to n. In particular, the value of Aj^n can be com- 
puted from the constant terms in the above formula. □ 

Corollary 2.11. For all n e , 

Jifn = ^ ] 7i,n/m 

for some real constants 7™^. 

Proof. Let i^' G be defined as in Lemma [2.91 Then JiKfn = Jif[n] is given 
by Proposition I2.10[ and the result follows from the commutation relations 
fl2.12ap - fl2.12ep and the invariance of the partial order under the Weyl group 
(cf. Lemma [231). D 

Proposition 2.12. For all n G there exist some real constants such 
that 

where 

En = J2 {'^Hi + b + b' + 2a{N - i)Y (2.18) 

i 

only depends on [n]. 

Proof. By Proposition 12.61 if = Ji is .^-invariant, so Hfn = K HKfn- By 
Lemma 12.91 one can assume that Kf^ = /[n] • By Proposition 12.101 

Hf^ = KHf[n^=Kj2-J!fln] 

i 

i j,m-<n j, m-<n 
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By Lemma 12.91 and Corollary 12.111 it suffices to show that ^ ^ is given 
by Eq. dnSD. 

There is no loss of generality in assuming that n G [Z^]. To compute 
the sum, let us gather those components of the multiindex that have the 
same value. Thus, take those components of the non-negative, nonincreasing 
multiindex n such that rik-i rik = ■ ■ ■ = rik+p ^ ^fc+p+i- Then i{nk+j) = k 
and ^{rik+j) = p + 1 for all ^ j ^ p. If rik ^ 0, 

Afc+j,n = 2nk+j+b+b'+2a{N-k-p+j) = 2nk+p-j+b+b'+2a{N-{k+p-j)) , 
so that 

k+p k+p 

Yl = Yl (2^^ + b + b' + 2a{N - i)f . (2. 19) 

i=k i=k 

If Uk = 0, Proposition 12.101 shows directly that Eq. (12.191) also holds in this 
case. □ 

Corollary 2.13. The spectrum of H as a self- adjoint operator on L^(T^) is 
given by {E^ : n G Z^}. Thus the degeneracy of each eigenvalue is at least 
\W\ = 2^N\. 

Proof. It follows from the density of the flag (I2.16P and the fact that H\-^k 
is a triangular matrix in the basis B with diagonal elements given by E^- □ 

Remark 2.14. Here we loosely (but safely) regard the spectra of the models 
H^^i as "sets with multiplicities". We shall use this convention whenever 
appropriate. 

Our goal now is to extend Corollary 12. 131 to the spin and spinless Hamil- 
tonians H^^i and ifsc- For simplicity, let us start with the latter, considered 
as an operator on L^(C). 

Theorem 2.15. spec(i/sc) = {^n : n G [Z^]}. 

Proof. By definition of the Weyl alcove, |2n|^^ times the inclusion C C 
yields a unitary transformation U from L'^{C) to the space of .^-invariant 
functions L'^{T^)^. Its inverse is given by |2n| times the symmetric extension 
iIj{w:k) = ip{'x) (for all w G W{BCn))- Since the action of and H agree 
on ^-invariant functions, H\i2^fNyr^ and H^c are isospectral. 

For each n G Z^, the orbit W{BCiy) ■ n intersects [Z^] exactly once, at 
[n] . By Lemma 12.91 this implies that 

{/„: |n| ^ A;, nG [Z^]} 
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is a basis of [Ti^)^, yielding a Hilbert basis of L^(T^)-'^ via the Gram-Schmidt 
procedure. By Proposition 12.121 the matrix of ifj^^fe^^ is triangular in an 
ordered basis and its eigenvalues are given by E^, where now n G [Z^]. □ 

The computation of the spectrum of H^^/ is similar, although the charac- 
terization of the bases is considerably more involved. 

Lemma 2.16. A basis of K^^i{l-L^ ® S) is given by 

= {A«'(/n|s)) : n G [Z\ |s) G S^'O ' (2-20) 



where 

Proof. Since 
it is clear that 



Si — Sj ^ when rij = Uj and i < j 

Si ^ |^-i,e' when rij = 



n^A,,.^ = X..'HA«'^, (2.21) 



{A,,,(/„|s)) : n G [Z% |s) G B^} 

must span the whole A^^'{T-l^ ® S). Moreover, a state of the form A^^K/nls)) 
(|s) G By,) vanishes if and only if some of the following conditions hold: 

(i) Hi = Uj, Si = Sj and e = —1. 

(ii) Hi = Uj = 0, Si = —Sj and e = e' = —1. 

(iii) Hi = 0, Si = and e' = —1. 

It is easy to check that the above states do not vanish and are in fact linearly 
independent modulo kerAg^'. Furthermore, it is immediate to see that for 
any element A^^/{fa\s)) not satisfying any of the above conditions one can 
use permutations and sign reversals to construct an element w G W{BCn) 
in the stabilizer of n such that n^[Aee'(/n|s))] belongs to the basis (12.201) . 
By Eq. (JMU), the claim follows. □ 

Corollary 2.17. For each element A^^/{fn\s)) of the basis B^^,, 



where 




if Hi > 
if Hi = . 



M+ = [MJ + 1 , M_ = \M] , (2.22) 

and [xj and [x] respectively denote the integer part of x and the smallest 
integer greater than or equal to x. 
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Theorem 2.18. spec(ifee') = {^n : n G [Z^], |s) G B^"'"']. 

Proof. By Proposition -f^ee' and H^'^' are isospectral. We can assume that 
the basis (12.201) of A^^/{T-C^ (g) E) is ordered according to the partial order -<. 
By Lemma [2.51 

//"Xe'(V^|s))=A«'((:ff^)|s)) 

for each G ^^(T^). Now the same arguments used in the proofs of Corol- 
lary 12.131 and Theorem 12.151 show that the action of H^^/ in the basis given 
by Lemma [2.161 is triangular with diagonal elements E,^. □ 



2.5 The spin chain 

We define the BCn spin chain Hamiltonian as 

8 H,,, = [ - Q (1 - eSi,) + sin-2(e, + Q (1 - e^)] 

+ 5^ (/3 sin-2 + /3' cos-2 eO (1 - e'S,) , (2.23) 

i 

where the spin sites 0<^i<---<^Ar<7r/2 satisfy 

- f3 cot ^i = Y, [cot(ei - e,) + cot(e, + Q] . (2.24) 

Remark 2.19. A cursory glance at the spin chain (11.270 reveals that the 
Hamiltonian (12.231) is in fact the appropriate BCn version of the celebrated 
HS chain. 

In this section we shall study the Hamiltonian (I2.23P and compute its 
spectrum by taking the large coupling constant limit of the model ( ll.lSdp . 
To this end, let us rescale the parameters of the BCn Sutherland Hamiltonian 
as 

b = aP, h' = aj3' . (2.25) 
The ground state function (12.31) can thus be written as 



where the smooth function 

A(x) = [logsin(xi — a;j) + logsin(xi + Xj)] [/? log sin Xj + log cos x^] 

i<j i 

(2.26) 
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is independent of the coupling parameter a. It is immediate to show that 
Eq. (12 .24^ merely states that ^ = . . . , ^n) G C is a critical point of A. In 
the following proposition we shall prove an existence and uniqueness result 
for the solutions of this equation. 

Proposition 2.20. The function A G L'^iC) fl C°°(C) has a unique critical 
point ^ in C , which is a hyperbolic maximum. 

Proof. Since C is bounded and A tends to — oo at dC, A must attain its 
maximum in C. Uniqueness shall follow easily from the fact that the Hessian 
D^X = (d'^X/dxidxj) is negative definite in the convex domain C. To show 
the definiteness of D^X, let us compute 



dx'j 



El 



sin ^(xj — Xj) + sin ^(xj 



+ Xj)~\ — /3 sin ^ Xj — P' cos ^ Xj < 



d^X 

dxjdxi 



— sin ^(xj — Xj) + sin ^(xj + x,- 



By Gerschgorin's theorem [96, Theorem 15.814], the eigenvalues of D'^X{x.) 
lie in 



u 



dxi 



with 



li 



dx\ 
d^X 



As d'l^X + 7i ^ — /3 sin ^Xj — /3'cos ^Xj < 0, D^A(x) is negative definite. 
Hence, any critical point of A must be a hyperbolic maximum. 

To complete the proof of uniqueness, let us suppose that and ^2 ^-^e two 
critical points of A. Consider the segment 7 : t G [0, 1] t-^ + (1 ^ ^ 

and set / = A o 7. Clearly /'(O) = /'(I) = 0, so there must exist t G (0, 1) 
such that 

r(t) = (7'(t),/^^A(7(t))y(t)) = o, 

contradicting the fact that D^X is negative definite in C. □ 

Remark 2.21. It can be proved [116,182] that ^ is also a (global) minimum 
of the part of the potential Vsc quadratic in a, namely 



C/(x) 



sm ^(x,— Xi 



)+sin ^(xj+Xj)+/5^ sin ^ Xj+/3'^ cos" 



Xj 
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Let us now consider the positive function 



J2a\ 



Ic 



^2aX 



e L^(C)nC°°(C). 



(2.27) 



c 



It is intuitively obvious that (po concentrates near its maximum ^ as a — > oo. 
More precise, one can prove the following 

Lemma 2.22. Let ^ be the unique maximum of X. Then (fo converges to 5^ 
in the sense of distributions and falls off exponentially fast away from ^. 



Proof. We shall show that 



lim / (/7o0 = 0(^) 
Jc 



for any e C^(C). By continuity, for any £ > there exists S > such that 

A(0-A(x)<£, |0(O-0(x)| <£, 

for all X in the ball B centered at ^ and of radius 6. Furthermore, one can 
assume that X\c\b ^ ~f with c ^ Cq > 0. Since the ground state function 
II is only defined up to a multiplicative constant and the claim obviously 
applies to constant functions, wc can assume that A(^) = 0(£) = and 
maxc \ 4>\ ^ 1. The method of steepest descent ensures that 



c 



TlV" 1 + 0(1) 



a J det[L>2^(^)] 



as a ^ oo, so that 



C\B 



Hence 



lim 

o— *oo 



e^"^ ^ Cia^/'e-"'^'^' = C»(a^/2g-ac52^ _ (2.28) 



C\B 



•^0 



c 



^ £ lim (po + lim / ^o\4>\ 
"■^°°Jb '^^°°Jc\b 

^ £ + lim C»(a^/2g-ac52^ 
a—*oo 



for any e > 0. 



□ 



Proposition 2.23. Let ip be a normalized eigenfunction of H^c- Then </? = 
IV^P — >• 5^ as a ^ oo and (p falls off' exponentially fast away from ^. 
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Proof. It is well known [104,109,159,160] that g = is a hypergeometric 
polynomial in z = (e^'^^ , . . . , e^'"^^) associated with the BCn root system. 
Therefore, is polynomially bounded in a and 

^ ^ 

C\B 

exponentially fast as a — cxd on account of Lemma [2.221 Since is obviously 
positive and J^ip = 1, the previous argument shows that 



lim / = 0(C) 



for all (j) e C^{C), as claimed. □ 

Remark 2.24. Finer control of the decay could have been obtained if needed [4, 
111]. 

We shall use the above lemma to relate the spectrum of the BCn chain 
with those of the BCn scalar and spin Sutherland models. The precise re- 
lationship can be introduced by means of the partition function of these 
systems, i.e., the trace of the exponential of minus their Hamiltonian oper- 
ators. In the large coupling constant limit, we obtain an exact formula for 
the eigenvalues of the BCn chain Hamiltonian. 

To this end, let us write the spin differential operator (ll.lSdl) as 

H,,: = ifsc + 8ahe,' , (2.29) 

where 

h£e'(x) = [sin~^(xj — Xj) (1 — eSij) + sin~^(xj + Xj){l — eSi, 



{p sin~2 Xi + P' cos~2 Xi) (1 - e'Si) (2.30) 



is an End(S)-valued multiplication operator and 

H.e' = h«'(l)- (2.31) 

Gloss 2.25. The essence of Polychronakos's freezing trick is that the eigen- 
functions of H^^/ become sharply peaked at ^ as a ^ oo, so that 

H,,,^{^) ~ i/sc^(x) + 8ah,,,(C)^(x) = (/fsc + 8aH,,0^(x) , 
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in view of Eqs. (12.291) and (12.311) . This idea suggests that the eigenvalue 
equation 

{H,,. - E,,)^< = 

should yield approximate eigenf unctions of H^c and H^e', i.e., 

(ifsc-^sc)^^O, (H,,, - E,,0^ - 

for some eigenvalues Esc-, E^e'. This leads to Polychronakos's heuristic for- 
mula [166] 

E,^, ~ Esc + 8a E,,/ 

relating the spectra of H^^i, H^c and H^^'- A major drawback of this relation 
is that it does not determine a priori which eigenvalues of the spin and scalar 
models should be combined to obtain an eigenvalue of the spin chain in the 
large coupling constant limit. This difficulty can be overcome [65, 84, 166] 
by summing over the spectra of these operators and promoting the above 
relation to the more precise formula (I2.34p relating the partition functions 
of i^ee', and H^^/ respectively. In what follows, we provide a full proof of 
this technique and apply it to compute the spectrum of the BCn chain. 

The partition function of a self-adjoint operator H acting on a Hilbert 
space Ti is the map Zh '■ K'*' — * IR'*' defined by 

Zh{T) = tin, e-^/^ = ^(^, e-^/^ij) (2.32) 

provided that the above sum converges (i.e., that e'^f^ is trace class for all 
T > 0.) Here B is an orthonormal basis of Ti. Clearly the above definition 
amounts to setting 

Zn{T)= J2 ^'^^^ 

Eespcc{H) 

when H has a pure point spectrum. We shall use the shorthand notation 
Z^^f = Zh^^,, Zsc = Zh^^ and Z^e' = Zy\^ , and define the equivalence relation 

A^ B lim — = 1. (2.33) 

a-^oo B 

Moreover, in the rest of this chapter we shall set E'o = by adding a constant 
term to the Hamiltonians H^^i and Hgc- The relation (I2.29P is not modified 
by this rescaling. 

Lemma 2.26 (Freezing trick). For all T > 0, 

lim Z^^i{8aT) 

Z^^,(T) = ^ . (2.34) 

lim Zsc{8aT) 

a— >oo 
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Proof. Let Bsc = {V^ = ^(x; a)} C L'^{C) and Sspm = {|s)} C S be orthonor- 
mal eigenbases of Hsc and H^e' respectively, and consider the orthonormal 
basis of H. given by the tensor product of these bases, i.e., B = = ip\s)}. 
As we have set i?o = 0, 

By Theorem 12. 18^ the eigenvalues of if^e' with £"0 = depend linearly on a, 
and hence it is not difficult to check that the sum defining ZH^^,/sa converges 
uniformly on compact sets. In particular, the limit 

lim Z^^/{8aT) 

a— >oo 

exists and depends continuously on T G M^. The same result holds for the 
scalar Hamiltonian H^c- 

Let us estimate the value of (\E', e~^"'/^""^\E'). By continuity and the fact 
that ^ is the unique maximum of A, for any e > there exist ei > and an 
open neighborhood i? of ^ such that 

-£< K,, - H,,/ < e (2.35) 

in B and A(^) — A > £1 in C\B. As H^^i ^ 0, it follows from Proposition 12.231 
that 



'2 vO 



Jc\B Jc\B 

as a —>■ 00. Therefore, 

lim(^,e-^«'/^"^^) = lim / e"^-'/^"^^) , (2.36) 

a— >oo a^oo J ^ 

and similarly for H^c + ca {c E M) since it only differs from H^^i by a times a 
relatively bounded operator. Moreover, Eq. (12.351) implies that 



Jb Jb Jb 

and hence 

lim (^, e-^"'/^'^^^) = (s|e-^^«'/^|s) lim (V^, e'^-Z^'^^V^) (2.37) 
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for each \I/ G i3, on account of Eq. (12.361) . As we have uniform convergence 
in compact sets we can exchange hmits, so that 



hm Z,,/(8aT) = hm e"^-'/^''^^) 

a— »oo a— ►oo ' ^ 

= y hm(^,e-^"'/«'^^^) 



y (s|e-^^«'/^|s) y hm(^,e-^-/8ar^) 
Z,,,{T) hm Zse(8aT), 



as claimed. 



□ 



We shaU now compute the partition functions of H^^i and ifsc and obtain 
the spectrum of the spin chain ( I2.23P using the freezing trick. A first obser- 
vation is that the eigenvalues 02.181) can be written in a more compact form 
as 

E„ = a^Eo + 8a ^ n,(p + N - i) + 0{l) , (2.38) 

i 

where n G [Z^], = \{(3 + (3') and Eq was defined in Eq. (El). 
Proposition 2.27. Set q = e~^/^. Then 



lim Zsc(8d 

a— >oo 



l-g 



Zsc(8aT) = 5^ e- 



ni(/3+Af-i) 



Proof. By Theorem 12.151 and Eq. (12.3 

neli'^] ne|: 
Defining Pi = rii — rij+i for 1 ^ z ^ — 1 and = nj^, we have 

J _ 



and hence 



n 



N i 



1-q 



ilf3+N-hi+l)] 



1 -1 



(2.39) 
□ 
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The computation of the partition function of the spin model is consider- 
ably more involved. First, let us define the set of partitions of a natural 
number n as 



^n = |J{keN'-:|k|=n}, 



r=l 



where |k| = fci + ■ ■ ■ + /c^ • We find it convenient to represent each positive 
nonincreasing multiindex n G [Z^] as 



n = (mi, . . . ,mi,m2, . . . ,m2, . . . ,mr, . . . ,mr), 



(2.40) 



where mi > m2 > ■ ■ • > m.,. ^ and ki = ^{nii). In particular, k G HW 
and condition (12.171) is satisfied if e = — 1. 

Example 2.28. = {(3), (2, 1), (1, 2), (1, 1, 1)}. Forn= (6,3,3,3,2,1,1) G 
[Z^ one has m = (6, 3, 2, 1) G [Z^] and k = (1, 3, 1, 2) G ^Pt n 

Proposition 2.29. Given k G ^n, set 



N, ^ iV,.(k) = ( 5^ fc,) + iV - - - - 5^ , (2.41) 

i=l i=l 

where 1 ^ j ^ r and r = r(k) is the length of h. Then one can write the 
large a limit of the partition function of the Hamiltonian (ll.lSdl) as 



lim Z^^SaT) = \ T\ 



2M+ l + (5i,,(% - 1)\ q 



kj 



1-g 



M,, + 5i^,{K - l)^ ^ /2M + 1 + 5i,,(A;, - \)\ g^-- 



1 - g^'- 



(2.42) 



Proof. It follows from Eq. fl2.40p that the eigenvalues (12.381) can be written 
as 

r feiH f-fci-i+fci 

^n~8aj]m, ^ (;^ + iV-j) 
i=\ j=fci+---+fe,_i+i 

r 



8a 5] m.fc//3 + iV - i - I - ] 
j=i V j=\ / 

r 



5 a mil/,; 
1=1 
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since we have set Eq = 0. The degeneracy (iee'(n) = dee'(k, m) of each eigen- 
value En is given by the cardinal of the set yB^'*^'*^ introduced in Lemma [2. 161 
It is not difficult to check that ^^'(k, m) solely depends on whether rrir = 
or not, and is in fact given by 



Ur^'^lY'^'^-'^) ^ dtAk) , ifm.>0, 
ci«'(k, m) = { ^"^ 

^M,+5,Akr-l)^^ j-j ^2Af+l+5,„(/c,-l)^) ^ ^o^^^j^) if = . 

The partition function can thus be written as 

nG[ZJV] 

r 



r-1 

rriiUi 



E [ E n + E (k) n ^ 

kG^Piv mi>--->mr>0 i=l mi>-->mr-l>0 i=l 



Since 



E n«-'"=En«--'"=Enn«''"- 

mi>--->ms>0 i=l peN« i=l pgN'' i=l i=j 

=nE^ 

pGN" j=l j=lpj>0 



we easily obtain Eq. (12.421) . □ 



From Propositions 12.271 and 12.291 we immediately derive the following 
theorem, which is the main result of this chapter. 
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Theorem 2.30. The partition function of the spin chain f l2.23p is given by 



N 



^AT) = n 



i=l 



l-q 



i[/3+JV-|(i+l)] 



ke<Piv 



Nr 



kr I [ kr Jl- g^'- 



r-1 

n 



2M + 1 + 5i,,ikj - 1)\ q 



kj J 1- q^^ 



. (2.44) 



Remark 2.31. The eigenvalues E of the spin chain Hamiltonian and their 
degeneracies d{E) = d^^'i^) can be recovered by identifying 

Z,AT)= Yl ^(E)g^ (2.45) 

EGspec(H^^/) 

where the sum runs over distinct eigenvalues. 



2.6 Statistical properties of the spectrum 



Several remarkable properties of the spectrum of the spin chain fl2.23p can be 
inferred from Theorem 12.301 First of all, for half-integer spin the partition 
function fl2.44p does not depend on e', since in this case M± = M + Hence 
the spectrum of the spin chain is independent of e' when M is a half-integer, a 
property that is not apparent from the expression of the Hamiltonian (12. 230 . 
Secondly, all the denominators 1 — q^^ , 1 ^ j ^ r, appearing in Eq. (I2.44p 
are included as factors in the product in the first line. Hence the partition 
function fl2.44p can be rewritten as 



where is given by 



<5e{o,i}^ 



^ - 1 
E, = ^^5,(/3 + iV--(z + l)), 

i=l 



(2.46) 



(2.47) 



and the degeneracy factor d^^>{M,S) is a polynomial of degree N in M. 
Therefore, 
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Proposition 2.32. For all values of (3, e, e' and M , the following statements 
hold: 

(i) For half -integer spin, H^e' and y\e-ti are isospectral. 

(ii) spec(Hee') C {E^ : d G {0,1}^}, and it exactly coincides with this set 
(of cardinal at most 2^ ) for sufficiently large M. 

(iii) Set Einax = I N{N + 1){2N + 3]3-2). Then 

H«' = Emax — ^-€,-€' (2.48) 

Proof. The first two statements follow from the above discussion. Let us 
prove the last claim. First, observe that the above equation holds with the 
constant 

4E^^ = J2 [sm-\^,~Q+sm-\^, + Q]+J2 (/? sin'^ + cos'^ e.) 

i^j i 

(2.49) 

by the definition of the spin chains. Since H^^i ^ and ker(Hee') = Aee'(S) 
is nontrivial for sufficiently high M, Emax must coincide with the highest 
possible eigenvalue of H_e which can be easily obtained using Eq. (12.471) : 

^ - 1 1 

Emax = ^(/^ + ^ - 2 + = 6 + ^^^^^ + ~ • 

i=l 

Note, however, that in the antisymmetric case some of the possible energies 
may not be attained if M is kept fixed and N increases. For instance, it is 
obvious that the kernel ker(H_^£) = A_ ,;(S) becomes trivial in this case. □ 

Remark 2.33. Eq. (I2.48P allows us to restrict our attention to the ferromag- 
netic chains H_^± without loss of generality. Moreover, E^ax = max spec (Hee') • 

We shall now present several concrete examples where we analyze the 
spectrum of the chains by means of Theorem 12.301 

Example 2.34. The structure of Eq. fl2.44p makes it straightforward to 
compute the spectrum of the spin chains for any fixed number of particles 
as a function of the spin. For instance, for = 3 sites and integer M the 

eigenvalues of the spin chain H are 0, /5 + 2, 2/? + 3, 3/? + 3, 3/5 + 5, 4/? + 

5, 5/3 + 6, 6/? + 8, with respective degeneracies 

|M(M- l)(M-2), |M(M2-1), |M(M + 1)(11M-2), 

|M(M+ l)(7M-4), I M(M + 1)(7M + 11), | M(M + 1)(11M + 13), 

I M(M + 1)(M + 2), |(M + 1)(M + 2)(M + 3). 
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Note that in this case spec (if ) = {E^ : <5 G {±l}}forM ^ 3, in agreement 

Proposition 12. 32[ The above results have been numerically checked for M = 1 
and a few values of (3 and jS' by representing the operators Sij and Si as 27x 27 
matrices. The obtained results are in complete agreement with those listed 
above. 

For a fixed value of the spin M, it is not apparent how to derive an 
explicit formula expressing the eigenvalues and their multiplicities in terms 
of the number of particles N. We shall next present two concrete examples 
for the cases M = 1/2 and M = 1. 

Example 2.35 (Spin 1/2). By Proposition 12. 32[ H_+ and H are isospec- 

tral. We have computed the partition function Z_ ± for up to 20 particles. 
For instance, for = 6 the antiferromagnetic spin chain eigenvalues and 
their corresponding multiplicities (denoted by subindices) are 

(9:^ + 32)2, (10;5 + 36)2, (li;5 + 38)2, (li;5 + 41)4, (12;5 + 38)i, (12:^ + 43)6, 
(13:^ + 43)3, (13:^+46)6, (14:^ + 46)4, (14:^ + 50)4, {15^ + 47)2, {15^ + 50)3, 
(16^+51)2, {10^ + 55)5, (17;5+53)i, (17;5 + 56)4, (18;^+ 58)3, 
(19^+61)2, (20;5 + 65)i, (2i;5 + 70)i. 

The number of energy levels increases rapidly with the number of particles N. 
For example, if = 10 the number of distinct eigenvalues (for generic values 
of /3) is 136, while for = 20 this number becomes 7756. It is therefore 
convenient to plot the eigenvalues E and their corresponding degeneracies 
d{E), as is done in Fig. l2.1l for = 10 particles. Note that Eq. f l2.47p implies 
that when /? ^ A^ the levels cluster around integer multiples of /3. In fact, 
for all A^ up to 20 we have observed that these integers take all values in a 
certain range jo, jo + 1, • • • , A^(A^ + l)/2. For example, in the case A^ = 6 
presented above, jo = 9. 

Example 2.36 (Spin 1). We have computed the partition functions Z_ ± of 
the spin chains H_^± with spin M = 1 for up to 15 particles. As remarked in 
the previous section, for integer M the partition functions Z_^± are expected 
to be essentially different. This is immediately apparent from Fig. 12.21 where 
the energy spectra of the even and odd spin chains Z_ ± with /? = \/2 for 
A^ = 10 particles are graphically compared. However, the standard deviation 
of the energy is exactly the same for both chains. This rather unexpected 
result will be relevant in the ensuing discussion of the level density (see 
Conjecture 12.391 below). We also note that, just as for spin 1/2, for A^ up to 
(at least) 15 and /? S> A^ the energy levels cluster around an equally spaced 
set of nonnegative integer multiples of p. 



200 



5 



10 



15 



20 



d{E) 



Figure 2.1: Eigenvalues E and degeneracies d{E) of the spin 1/2 
chain H_ -t for = 10 particles and P = \/2- 



The numerical explorations that we have carried off naturally give rise to 
several conjectures that we shall now discuss in detail. 

Conjecture 2.37. For f3 ^ N, the eigenvalues cluster around an equally 
spaced set of levels of the form i(3, with j = jo, jo + 1, . . . , A^(A^ + l)/2 

Remark 2.38. For sufficiently large values of the spin M this conjecture (with 
jo = 0) follows directly from Eq. fl2.47p . Numerical calculations for a wide 
range of values of and M fully corroborate the above conjecture. 

Conjecture 2.39. For iV 1, the level density follows a Gaussian distri- 
bution 

More precisely, we claim that the number of eigenvalues (counting their 
degeneracies) in an interval / C M is approximately given by 




(2.50) 



where 
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Figure 2.2: Comparison of the energy levels E and degeneracies 

(i(E) of the spin 1 chains H (left) and H_+ (right) for = 10 

particles and (3 = 



is the normal (Gaussian) distribution with parameters and a respectively 
equal to the mean and standard deviation of the spectrum of the spin chain. 
Although the shape of the plots in Figs. 12. I land l2.2l make this conjecture quite 
plausible, for its precise numerical verification it is preferable to compare the 
distribution function 




(2.52) 



of the Gaussian probabihty density with its discrete analog 

F(E) = (2M + 1)'^ ^ d(e), (2.53) 

where (i(e) = (i«'(e) denotes the degeneracy of the eigenvalue e. Indeed, our 
computations for a wide range of values of M and > 10 are in total agree- 
ment with the latter conjecture for all four chains (12.231) . This is apparent, 
for instance, in the case /? = v/2, M = 1/2, and A^ = 10 presented in Fig.EJl 
The agreement between the distribution functions (12.521) and (12.531) improves 
dramatically as A^ increases. In fact, their plots are virtually indistinguish- 
able for A^ > 15. 
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Figure 2.3: Distribution functions Fj^{E) (continuous line) and 
F(E) (at its discontinuity points) for (3 = \/2, M = 1/2, and 
= 10. 



It is well known in this respect that a Gaussian level density is a char- 
acteristic feature of the embedded Gaussian ensemble (EGOE) in Random 
Matrix Theory [152]. It should be noted, however, that the EGOE applies 
to a system of particles with up to n-body interactions {n < N) in the 
high dilution regime N ^ oo, k — oo and N/k, 0, where k is the number 
of one-particle states. Since in our case k = 2M -|- 1 is fixed, the fact that 
the level density is Gaussian does not follow from the above general result. 

If Conjecture 12.391 is true, the spectrum for large N is completely char- 
acterized by the parameters fi and a through the Gaussian law (12.511) . It 
is therefore of great interest to compute these parameters in closed form as 
functions of and M. To this end, let us write 

where the constants h^, and hi can be easily read off from Eq. (12.231) . 
We shall denote by fi± the average energy of H_ j-. 
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Table 2.1: Traces of products of the spin operators. 



Operator 


Trace (inte^ 


^er 


M) 


Trace (half-integer M) 


S^ 


(2M+ 1)^- 


-1 







Sij ; Sij 


(2M+ 1)^- 


-1 




(2M+ 1)^-1 


SiSj 


(2M+ 1)^- 


-24 


-2<5ij 


(2M + l)^6ij 


SijSk, SijSk 


(2M+ 1)^- 


-2 







Sij Ski 


(2M+ 1)^- 


-2 






Sij Ski, Sij Ski 


(2M+ 1)^- 


-24 


-2<5ifc<5ji-l-25ii5jfc 


(2M + l)^-2+2<5ifc5ji+2<5,j5jfe 



Proposition 2.40. // i/ie spin M is an integer, 



= 6(2m"+\) + ^^^^^ + 3/5 - 2) , (2.55) 



1 



2M + 

wi/i Si = 'Ylii^i- If it is a half-integer, 

1 



- [(M + l)E^ax-2Si] , (2.56) 



2M + 



- [(M+l)E^ax-Si] . (2.57) 



Proof. We shall begin with the case of integer spin. Using the formulas for 
the traces of the spin operators given in Table 12.11 we immediately obtain 

/i_ = (2M + l)-^trH__ 

2(M + 1) r^,^ ~ >^ 1 M+1 ^ 

i^j I 



6(2M + 1) 

Here we have used Proposition 12.321 and Eq. fl2.49p to compute the sum. On 
the other hand, the average energy /x+ of the chain H_+ is given by 

^+ = (2M + l)-^trH_+ 

2(M+ 1) ~ 2M 
= ^ > {hii + hii) H > hi 

i^j I 
= [(^ + l)Emax - 2Si] , M integer, 
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where Si = ^ . hi. 

When M is a half-integer, /i+ = /i_ due to Proposition 12.321 A similar 
calculation shows that the formulas for the traces of the spin operators in 
Table 12.11 yield the expression (12.571) for the mean energy. □ 

Let us turn now to the (squared) standard deviation of the spectrum of 
H_ ±, which is given by 

2 tT{Hl^^) (trH_,±)2 



^ (2M + 1)^ (2M+1)2^' 

A long but straightforward calculation using the formulas in Table fTTl yields 
the following 

Proposition 2.41. For integer spin, the standard deviation is given by 

4M(M + 1) 

^± - (2M + 1)2 ' ^2-^^^ 
with S2 = 2 '^{hfj + h'^j) + ^hf. If M is a half-integer, then 



4M(M + 1) 
(2M + 1)2 ' M(M + 1) 



(2.59) 



where ^3 = i Yl - hijhij. 



Since Si, S2 and S3 do not depend on M, the above propositions com- 
pletely determine the dependence of /i± and cr-t on the spin. A nontrivial 
consequence of Proposition 12.411 is the equality of the standard deviation of 
the energy for the even and odd antiferromagnetic chains (for half-integer 
spin, this trivially follows from the fact that the even and odd chains have 
the same spectrum). This result is quite surprising, since for integer spin 
the energy spectra of the chains H_^-|- are essentially different, as shown in 
Fig. O 

To analyze the dependence on the number of particles of the spectrum 
we still need to evaluate Si, S2 and S3 as functions of A^. Although it is not 
at all apparent how to compute these quantities in closed form, in view of 
Eqs. (I2.55p -f l2.57p one is lead to formulate a natural 



Conjecture 2.42. The average energy n± and its squared standard deviation 
a\ depend polynomially on N. 
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In fact, since Emax is a polynomial of degree 3 in iV by Proposition 12.321 
it follows that the degrees in N of /i± and cx^ cannot exceed 3 and 6, respec- 
tively. The latter conjecture and this fact allow us to determine the quantities 
Ej(A^) by evaluating ^± and cr^ for = 2, . . . , 8 and M = 1/2, 1 using The- 
orem [2]30] (cf. Eqs. f l2.56p . (12.581) and (12.591) ). Hence one immediately proves 
the following 



Proposition 2.43. // Conjecture 2.42 is true, then 



51 = 8E3 = ^(2/3 + iV-l), 

/V _2 _ 

52 = —[2{2N^ + 3N + 13)/3 + {N - 1)(5N'^ + 7N + 20) p 

+ -{N-1){8N^ + 3N^ + 13N-12)]. (2.60) 
5 

These expressions, together with Eqs. (I2.55p - (l2.59p . completely determine 
/i± and cr± for all values of M and A^. It has been numerically verified that 
the resulting formulas yield the exact values of /x± and a± computed from 
the partition function (12.441) for a wide range of values of M and N. This 
provides a very solid confirmation of Conjecture 12.421 Let us mention, in 
closing, that formulas analogous to (12. 550 - 02.591) expressing the mean and 
standard deviations of the energy for the ferromagnetic chains H_|_ ± can be 
immediately deduced from the previous expressions and Eq. (12. 48 p . 



2.7 Models with constant magnetic fields 

In this section we shall briefly show how a constant magnetic field can be 
added into the picture that we have presented in the previous sections. For 
the sake of simplicity, we restrict ourselves to the case M = 1/2, e = — 1 
("one-dimensional electrons"), and define the magnetic BCn Hamilto- 
nian as 

^mag _ ^ ^ _ g^^ _ ^2.61) 

i 

Gloss 2.44. The motion of a (three-dimensional) electron in the presence of 
static magnetic and electric fields B = V A A and E = — in is 
governed by the Pauli operator 

H={iV + eAf - 2eB ■ S + el/ , 

which is the nonrelativistic limit (in the resolvent norm topology) of the Dirac 
operator [193]. Here x = (x, ?/,2;), e is the electric charge of the electron. 
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S = |cr is the spin operator of the particle and cr = (cxi, (T2, as) are the Pauh 
matrices. If B is constant, one can take A = |x A B, so that 

e2 

if = -A + ieB ■ X A V + — llx A Blp - 2eB • S + el/ . 

4 II II 

Since we are interested in a one-dimensional scenario, we drop the second 
and third summands in order to preserve the symmetries dy and dz- This is 
a priori a good approximation to the problem when the angular momentum 
and the charge is small and the electric potential is strongly binding. 

Hence one is led to define the magnetic BCn Hamiltonian for spin 1/2 
particles as 

i 

As the exchange operator Sij is written in terms of the SU(2)-invariant prod- 
uct as in Eq. f ll.26p . in this model of one- dimensional particles there is only 
one privileged direction e, namely, the one which implements the Weyl re- 
flections in spin space via 

Si —— 2g ■ Sj • 

Therefore it is natural to choose B parallel to this direction, so that 

2e^B-Si = e||B||^5i. 

i i 

If we now rescale e||B|| to 8aB as in Eqs. (12.251) we arrive at Eq. (12.611) . 

Since if^e' and Si commute, they possess a common basis of eigen- 
functions. Let n G [Z^] and |s) G -B"'"'*^ (cf. Eq. Lemma [2.161) . and define 

p(k) = card {z : fcj = l} , 
p+(n, s) = card {i : St = |, rii ^ 0, #(ni) = l} , 

where k is defined as in Eq. f l2.40p . Then one can easily prove the following 
Theorem 2.45. Let B = B^''^ . Then 

spec(ifr') = {^n+8a5[p(k)-2p+(n,s)-5„^,o(l+e)] : n G [Z^], |s) G B] . 

Proof. Let Si = C\j) Q) C| — |) be the one-particle spin space, so that S = 
Sf^, and define the basis of one-particle states 

|±) = |i)±|-i)GSi. 
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Clearly 5*1 1±) = ±|±). Now define the one- and two-particle operators 

A«:|±|)^±|e), 
Af):|±|,Ti)-|e,-6), 

with e = ±. By Propositions 12.31 one can consider the action of the differen- 
tial operator flMTD on A_,,(L2(T^) ® S). Theorem [2IIH] states that H_^, is 
diagonal in the basis (12.201) . Expressing the nonnegative, nonincreasing mul- 
tiindex n in terms of m and k as in Eq. (12.401) . one can use Conditions (i)-(iii) 
in Lemma 12.161 to readily verify that 

(A_,,(/„|s)):|s)GS> = (^„,s:|s)Gi3>. 

Here we have set 



with 



and 



^n,s = A_,, ( /„ (g) v,(sk^+i, kj, rrij) j 



1=1 



v,(±i, l,m 



±) if m 7^ , 

e) if m = - 

v,(±i,2,m) = 1+)®!-), 



so that for each j = 1,2 one has 

j 

^S'iV,(±i,j,m) 

i=l 



±(5j^i if m ^ . 

e if m = - 



Since Yli^i commutes with the symmetrizer and rrij = n^^+i, one immedi- 
ately derives from these equations that 



J^<S'i^n,s = ^n,s 7] card {i : Si = ^, 0, = l] + e5, 



njv,0 



ri=±l 



= [2p+(n, s) - p(k) + 5„^,o(l + e)] ^n,s , 

i.e., Yli^i is diagonal in this basis. As the action of H_ ^ and H agree on 
\E'n s and, on account of Proposition 12. 121 the action of the latter operator on 
{/n} is triangular with diagonal elements E^, the result follows. □ 
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The spin chain associated with the magnetic Hamiltonian (12.611) is 

H^"s = H_,, - 5 ^ . (2.62) 

i 

It is straightforward to show that Lemma [226] also apphes, mutatis mutandis, 
to the magnetic Hamiltonian fl2.6ip . Hence we shall compute its spectrum by 
means of the freezing trick. To this end, let us denote by Z™*^^ and Z™^^ the 
partition functions of (12.611) and (I2.62p . and again we shall set the ground 
state energy Eq equal to 0. We shall freely use the notation of Section [231 

Proposition 2.46. Let 

N N 

«7V = U {k e {1, 2? : |k| = iV} = n |J{1, 2y 
the set of partitions of N by {1,2}. Then 

ke2Sjv i=i ^ ^ ^ ^ H / 

(2. 

Proof. From Proposition 12.291 one immediately obtains 

Zf'''^{8aT) ~ ^ ^ gELi'"^!'- ^ ^B[2p+(n,s)-p(k)] ^ 

kg23jv »ni>--->r?ir^O |s)gB 

so that the sum over |s) depends on whether m^. > or not. If this is the 
case, 

card{|s) e B : p+(n,s) = g} = (^^^^^^ (2-64) 
for any ^ ^) ^ p(k), so 

^i?[2p+(n,s)-p(k)] ^ ^-i?p(k) J2 fP^^Aq-^Be ^ (^^B ^ q^Byi^) _ 
|s)eB £,=0 ^ ^ ^ 

If, on the contrary, m,, = 0, then kj. = 1 and it is not difficult to see that the 
formula (I2.64p becomes 

card (Is) G B : p+(n,s) = g} = (^^^^|^" , 
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where ^ ^ p(k) — 1. Hence in this case 

p(k)-i 



^ ^iJ[2p+(n,s)-p(k)+l] ^ ^-i?[p(k)-l] ^ /^P(k) ^\q2Be ^ ^ ^_B^p(k)-l 



|s)eB e=o 
and therefore 

kG2Siv m.i>->m, — 1>0 

Taking Eqs. fl2.41l) and 02.431) into account we obtain the desired formula. □ 



The freezing trick formula fl2.34l) and the scalar partition function (cf. 
Proposition 12.271) then yield the following 

Theorem 2.47. When M = 1/2, the partition function of the spin chain (12.621) 
is given by 

2^^^(T) = JJ [1 - 5 (*+!)]] 

keosjv 



Chapter 3 



Spin models with NN 
interactions 

3.1 Introduction 

In this chapter we analyze in detail three QES iV-body spin models which 
can be regarded as cyclic versions of (generalized) CS models, as constructed 
by Finkel et al. [81]. As this kind of systems feature only near-neighbors 
interactions, we shall generically call them A^'A'' models. There are two rea- 
sons that make this kind of systems very promising from a physical point of 
view. First, some of them are related to the short-range Dyson model in ran- 
dom matrix theory [28] . Second, the HS chains associated with these models 
occupy an interesting intermediate position between the Heisenberg chain 
(short-range, position-independent interactions) and the usual HS chains 
(long-range, position-dependent interactions). We shall consider this latter 
aspect in detail in the next chapter. 

The study of these systems was initiated a few years ago, when Jain and 
Khare presented a novel class of partially solvable models in which each par- 
ticle interacted only with its nearest and next-to- nearest neighbors [124]. In 
a subsequent paper [10], Auberson, Jain and Khare discussed a generaliza- 
tion of these models to the BC^ root system and to higher dimensions. The 
latter papers, however, left open some important issues, such as the exact 
or quasi-exact solvability of these models, the derivation of general explicit 
formulas for their eigenf unctions, or even the existence of similar models for 
particles with spin. The last question was first addressed by Deguchi and 
Ghosh [48], who introduced spin 1/2 extensions of the scalar models of Jain 
and Khare and used the supersymmetric approach to obtain some eigenfunc- 
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tions. All these authors solely managed to construct a few exact solutions, 
with trivial spin dependence. Moreover, their approach was by no means 
systematic. 

Our construction of NN models is based on a nontrivial modification of 
the Dunkl operator technique and yields fully explicit formulas for a wide 
variety of spin eigenf unctions. The presentation of these results is organized 
as follows. In Section 13.21 we define the Hamiltonians of the spin iV-body 
models which are the subject of this chapter (cf. Eqs. fl3.ip -( l3^ ). and show 
that they can be expressed in terms of suitable differential operators with 
near-neighbors exchange terms. Section 13.31 is entirely devoted to the char- 
acterization of certain finite-dimensional spaces of polynomial spin functions 
invariant under the latter operators. In Section 13.41 we show that the eigen- 
value problems for the Hamiltonians of the models (13. 2p restricted to their 
invariant spaces reduce to finding the polynomial solutions of a correspond- 
ing system of partial differential equations. By completely solving the latter 
problem, we obtain several (infinite) families of eigenf unctions of the mod- 
els (13. 2p in closed form. 

The material presented in this chapter is taken from Refs. [66-68]. 



3.2 The NN models 

In this chapter we shall study in detail the A^-body NN spin models defined 
by the Schrodinger operators 

if, = -A + K, (3.1) 

where e = 0, 1, 2 and 

2a^ sr-^ 2a(a — Sj^j+i) 



Vo = u} r + > + > — , 3.2a 

^ {Xi - Xi-i){Xi - Xi+i) ^ [Xi - Xi+i)^ 

^ x? ^ (x2-xti)(a;2-xf+i) 



+ 4« E 7^^^(« - ' (3.2b) 



V2 = 1c? cot(xi - Xi_i) cot(xi - Xj+i) + -r-TF — '^'''"^^l , (3.2c) 
^ ^ sm [Xi - Xi+i) 

where we identify xq = xat and xtv+i = Xi. As we did in the study of CS 
models, we shall assume a, 6 > ^ and take into account the inverse-square 
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singularities to consider the above Hamiltonians as self-adjoint operators on 

where the "reduced" configuration spaces are given by 

Co = {x : Xi < • ■ ■ < Xat} , (3.3a) 
Ci = {x : < a;i < ■ ■■ < xjv} , (3.3b) 
Ca = {x : < Xi+i - Xi < IT, \fi < N] . (3.3c) 

Remark 3.1. Unlike CS models, the Hamiltonians defined by (13.21) on Ti.^ and 
on L^(R^) (8) S (or L^(T^) S for the third potential) are a priori different. 
This is due to the fact that the Weyl symmetry present in CS models has 
been replaced in the NN case by the cyclic group 

3 = (ni2---n^i> (3.4) 

~ (o"ei— 62 ■ ■ ■ ^ejv— ei) ) 

and this latter group does not tessellate (or T^). Nonetheless, all the 
eigenfunctions presented in Theorems 13.171 13.191 and 13.231 can be clearly ex- 
tended to eigenfunctions of H^ on the "unreduced" Hilbert space since all 
our developments are essentially algebraic. 

We shall also make use of the scalar reductions of these models, which 
are the self-adjoint operators on L^{C^) given by 

H^ = H,\s.,.^,^i ■ (3.5) 

The identity 

i/,(^is)) = (/7:»is), 

with |s) a symmetric spin vector, is obvious. 

The models (13. ip share a common property that is ultimately responsi- 
ble for their partial solvability, namely that each Hamiltonian H^ is related 
to a scalar differential- difference operator involving near-neighbors exchange 
operators. In fact, consider the operators 

T, = J2 4dl+'^ci J2 ^— i4d^ 

^ ^ Zi- Zi+i 

where di = S^., z^+i = Zi, and 



-Zi^-^d,+i)-2a2^—— -(l-Ki,i+i), 

^ {Zi Zi^i) 

(3.6) 



^?o(a;,2/) = 1 , ^i{x,y) = -{x + y) , ^2{x,y) = xy. 
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Table 3.1: Parameters, gauge factor and change of variable. 





e — u 


e — i 


e — z 


c 


-1 


-4 


4 


c_ 





-2(26 + 1) 





Co 


2u 


4cj 


4(1 - 2a) 


Eq 


Nu{2a + 1) 


Ara;(4a + 26 + 1) 


2Ara2 


/i(x) 


i 


_^ ^2 -p-|- 12 9 Ifi f) 

i 


Hsin" \xi - Xi+i\ 

i 




X 







Each Hamiltonian is related to a linear combination 

= cT, + c_J^ + CoJ° + ^o (3.7) 
of its corresponding operator T^ and the first-order differential operators 

J- = Y,^^. J° = 5^^,a, (3.8) 

i i 

through the star mapping (11.221) . a change of variables and a gauge transfor- 
mation. More precisely, 

H, = fiH:\^^^^^^^^fi-\ 6 = 0,1,2, (3.9) 

where the constants c, c_, cq, Eq, the gauge factor /i, and the change of 
variables ( for each model are listed in Table 13. 1[ Here and in what follows 
we write the star mapping and the symmetrizer as * = ^ and A = A+, since 
for the sake of concreteness we shall restrict our attention to the symmetric 
case. 

From Lemma [23] and Eq. flHUD it follows that if <l>(z) G A[C°° (g) S) is a 
symmetric (formal) eigenfunction of H^, then 

vl>(x)=^(x)$(z)|,,=C(x,) (3.10) 

is a (formal) eigenfunction of with the same eigenvalue. In the next section 
we shall construct a flag C Til C ■ ■ • of finite-dimensional subspaces of 
A(C[z] (§5 S) invariant under each H^. We will show that the problem of 
diagonalizing in each subspace is equivalent to the computation of 
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the polynomial solutions of a system of linear differential equations. We shall 
completely solve this problem, thereby obtaining several infinite families of 
eigenfunctions of H^ for each e. From the expressions for the change of 
variable and the gauge factor in Table 13.1^ and the fact that the functions 
$ in Eq. f l3.10p are in all cases polynomials, it immediately follows that the 
eigenfunctions thus obtained are in fact normalizable. 

Gloss 3.2. The above procedure parallels the construction of (generalized) 
CS models of A^ type [81], which also applies to the BCn case [82]. Each 
operator T^ is the NN analog of the sum of squares of one of the three known 
families of An Dunkl operators. E.g., the operator Tq is to be compared with 
the sum of the squares of the Dunkl operators (11. lip associated with the Ajy 
system, i.e.. 

Hence, one can regard Hq as the proper NN analog of the Calogero model 
from the point of view of quasi-exact solvability. In contrast with the usual 
CS models, however, the physical Hamiltonians H^ associated with the cyclic 
differential-difference operators T^ include a three-body term which turns out 
to be crucial in the analysis of the solvability properties of these models. 

Two natural questions arise from the above discussion. First, it should 
be ascertained whether the full Dunkl operator approach can be extended 
to the NN setting or one is actually obliged to do without it. Secondly, one 
could wonder whether one would be able to use the above construction to 
prove the quasi-exact solvability of the models (13. ip as in the case of CS 
models. The answer to both questions is negative. Since 3 is not a Coxeter 
group, there are no known procedures to construct a set of Dunkl operators 
associated with this problem. Moreover, a couple of technical details prevents 
the usual approach to CS systems to work successfully in the NN context. 
Let us elaborate on this. 

The proof of the quasi-exact solvability of the Hamiltonian obtained from 
Tq^ (or a linear combination similar to (13. 9p ) through a gauge transformation 
is based on two facts: 

(i) Tq^ leaves invariant the space of polynomials 

P"={/GC[z]:deg/^n} 

for any n G Nq, where here and in what follows we denote by deg / the 
degree in z of the polynomial /. 



64 



CHAPTER 3. SPIN MODELS WITH NN INTERACTIONS 



(ii) Tq^ commutes with the symmetrizer A, which satisfies KijA = SijA. 

These conditions immediately imply that 

(ToCS)*A(p"|s)) = A[{T^^V'')\s)] c A(P"|s)) 

for any n G No, |s) G S. Whereas Tq certainly satisfies the first condition, it 
does not commute with A, as it is not invariant under arbitrary permutations 
of the particles. It is also clear that any attempt of replacing A by the 
symmetrizer under the cyclic group (13.41) is doomed to fail, as this symmetry 
is not enough to exchange -R'i^j+i by ^j^j+i, and thus by H^. 

Hence, the extension of the CS methods to the context of NN models is 
not at all trivial, and in fact is not a priori guaranteed to admit finite- 
dimensional invariant subspaces of /iA(C[z] ® S). We shall next show that 
this is the case nonetheless. 



3.3 The invariant spaces 

In this section we shall prove that each operator leaves invariant a fiag 
C C • • • , where T^" is a finite-dimensional subspace of A(P" ® S). 

This result will then be used to construct a corresponding invariant fiag 
C Til C ■ ■ • for the operator H^, where C 7^" for all n. 

Let us first introduce the following two sets of elementary symmetric 
polynomials: 

ak = ^Zi, Tk= Zi,---Zi^; k = l,...,N. 

i ii<---<ik 

It is well known [145] that any symmetric polynomial in z can be expressed 
as a polynomial in either cr = (cri, . . . , ctat) or t = (ti, . . . , r^). 

We shall denote by 2aX^ the terms of T^ linear in derivatives, that is 

^e = y] {Z'idi - ^i+i^i+l) . 

Zi — Zi^i 

In the next lemma we show that each vector field leaves invariant a 
corresponding fiag C C • ■ ■ of finite-dimensional subspaces of the 
space A(C[z]) = C[cr] = C[t] of symmetric polynomials in z. If / is a 
function of the symmetric variables cxi, cr2, as, Tat-i, Tat, we shall use from 
now on the convenient notation 




d^J-, fc = l,2,3, 
drJ, k = N-l,N. 
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Lemma 3.3. For each n = 0, 1, ... , the operator leaves invariant the 
linear space X^, where 

= C[ai, ^2, a3]nP", = C[ai, ^2, r^]^P^ X^ = C[ai, r^_i, r^]nP". 

Proof. Let us first consider the vector field Xq. Since 

V Zi — Zj+i -^-^ — / 

if / G A'J* we have 

Xo/ = 2(iV/2 + 3(Ti/3) G • (3.11a) 
The proof for the remaining two cases follows from the analogous formulas 

Xi/ = iV/i+4ori/2, feX^- (3.11b) 

X2f = 2a Ji + N{rN-ifN-i + tn/n) , f e X^ . (3.11c) 

□ 

Remark 3.4. It should be noted that these fiags cannot be trivially enlarged, 
since, e.g.. 



2(T2 - 


1" ^ ] ^i^^i+l , 
j 






2(T2- 


i 


XiTn^I 


i 


2(T2- 


1" ^ ] ) 




= ntn-2 - tn y^^{ziZi+iy^ 



are not symmetric polynomials. 

We note that the restriction of T^ to X^ C A(C[z]) obviously satisfies 

T,\;,n = J2z'id^ + 2aX,. (3.12) 

i 

The second-order terms of the operator fl3.12p . however, do not preserve 
the corresponding space X^, unless one imposes the additional restrictions 
specified in the following proposition: 
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Proposition 3.5. For each n G No, the operator T^ leaves invariant the 
linear space 5", where 

= {/ G : /33 = 0} , 
5r = {/ G A-i" : /22 = fNN = 0} , 
^2 = {f ^ '^2 '■ /ll = fN~l,N~l = 0} . 

Proof. Let us begin with tlie operator Tq. If / G Xq, an elementary compu- 
tation shows that 

dJ = h + 2zj2 + Szfh (3.13) 

and therefore 

= ^(/ii + 2/2) + 2(2/12 + 3/3)^1 

i 

+ 2(3/i3 + 2/22)^x2 + 12/23^3 + 9/33^4 . (3.14) 

From the previous formula and Eq. (13.1 lap it follows that Tq/ G Sq whenever 
/ G 5q . Similarly, if / G Af" we have 

dif = A + 2zj2 + z-^^fM , (3.15) 

so that 

^^^'f = (/ll + 2/2)0^1 + 4/120^2 + 4/22^3 

i 

+ 2NfiNTN + 4/2ArairAr + /nnTn-iTn , (3. 16) 

which together with Eq. (iS.llbll implies that Ti/ G 5" for all / G 5". Finally, 
if / G then 

dif = fl+ (^r^^Af-l - Zl'^TN)fN-l + Zi^TNfN (3-17) 

and hence 

Y^l^U = /ii^2 + 2/i,^_i((Tirjv-i - Ntn) 

i 

+ 2/iArCrirAr + f^-l^N-l [{N - l)r^_i - 2r7V-2TAr] 

+ 2(iV - l)fN-l,NTN-lTN + iV/iV7v4 • (3-18) 

The statement follows again from the previous equation and Eq. fl3.11cp . □ 
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The last proposition implies that each operator preserves "trivial" 
symmetric spaces 5"® A(S) spanned by factorized states. The main theorem 
of this section shows that in fact the latter operator leaves invariant a flag 
of nontrivial finite-dimensional subspaces of A{V"' ® S). Before stating this 
theorem we need to make a few preliminary definitions. Given a spin state 
|s) G S, we set 

I^^) = mE^I^)' ''^^■^ = M S ^(1±5i2)|s), (3.19) 

7r(l)=i 7r(l)=j,7r(2)=j 

where Sn is the symmetric group on elements and we identify an abstract 
permutation vr with its realization as a permutation of the particles' spins. 
From Eq. (13.191) we have 

A(/(2i)|s)) =^/(zO|s«), A{g^iz,, Z2)\s)) = ^^±(z„ z,)!^) ' 

i i<j 

(3.20) 

where the last identity holds if g±{x,y) = ±g±{y,x). We also define the 
subspace 

S' = { |s) G S : EKm) e A(S) I C S . (3.21) 
Theorem 3.6. Let 

V = {f{(^i,(^2,CT3)A\s),g{ai,a2,CT3)A{zi\s)),h{ai,a2)A{zl\s)), 

h{ai,a2)A{ziZ2\s')),w{cri,a2)A{ziZ2{zi - Z2)\s)) : /sa = 5(33 = O) , 
^1" = (/(o^i,o^2,TAr)A|s),5(((Ji,rAr)A(2;i|s)) : f22 = fNN = 9nn = 0) , 

V = (/(cri,T-7V-l,rAr)A|s),5((rAr_i,rjv)A(zi|s)), 

TNq{cri,TN)A{z^^\s)) : fu = fN-i,N~i = gN-i,N-i = qu = 0) , 

where |s) G S, |s') G S', deg/ ^ n, deg^f ^ n — 1, degh ^ n — 2, deg/i ^ 
n — 2, deg u; ^ n — 3, deg g ^ — + 1. T/ien TJ^ is invariant under T^ for 
all n G No. 

Proof. By Proposition 13.51 it suffices to show that T^ maps / (iS" ® A(E)) 
into T^. We shall first deal with the operator Tq. Consider the states of the 
form gA{zi\s)), with g G ^o""*^. Since 

{di - di+i)zi = ^ (1 - Ki^i+i)zi , V i, /, 
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we have 

Toigzi) = {Tog)zi + 2dig . 

Calling 

$(fc) ^ $W(|s)) = A{z'[\s)), keZ, (3.22) 
from Eqs. fl3.13p and (13.201) we obtain 

3 

M9^^'^) = J2'^o{9z^)h) = (To(7)$« + 2 5^A;^^,$(^-i) G . (3.23) 

i k=l 

Similarly, if h{ai,a2) E Sq'"^, the identity 

{di - di+i)Zi = ^ (1 - Ki^i+i)Zi + {zi - zi+i){5ii + 5/,i_i) , Vi, / 

zi - Zi+i 

implies that 

To{hz}) = {Toh)z^ + AzAh + 2(2a + l)h , 

and therefore 

To{h^^'^)=J2Toihz^)h) 

i 

= {Toh + 8/12)$^^^ + 4/ii$(^) + 2(2a + l)/i$(°) (3.24) 

belongs to '7^"~^ on account of Eqs. (I3.11al) and (13.141) . On the other hand, 
from the equality 

{di - di+i)ziZj = (1 - Ki^i+i)ziZj - {zi - 2;;+i)5j.i+i5i,i, Mi <i,\ll 

Zl — Zij^x 

it follows that 

Ti^ihziZj) = (TQh)ziZj + 2{zidjh + Zjdih) — 2ah6j^i+i . 

Setting 

$(2) = $(2)(|s)) = A(^i^2|s)) (3.25) 
and using again Eqs. (13.131) and (I3.20p we then have 

To{h¥'^) =J2ToChz.z,)\sp 

i<j 

= {Toll + 8/12)$^'^ + 2hiA[{zi + Z2)\s)] - 2ahJ2Ki+i) ■ (3-26) 
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Since A[{zi + Z2)\s)] = A[zi{l + 5i2)|s)], the RHS of Eq. (Km belongs to 
7^"~^ if and only if |s) G S'. The last type of states generating the module 
7q" are of the form w{ai, a2)^^^\ where 

$(3) ^ $(3)(|s)) ^ AizM^i - ^2)|s)) . (3.27) 

From the equality 

1 1 

{di - di+i) [ziZj{zi - Zj)] = ■- (1 - Ki^i+i) [ziZj{zi - Zj)] 

Zl — Zi+i [Zi — Zi^i) 

+ {Si^i-i + Sii)zj - {6ij-i + 6ij)zi, Wi < 

it follows that 

To[wZiZj{zi - Zj)] = {Tow)ziZj{zi - zj) + 2zj{2zi - Zj)diW 

- 2zi{2zj - Zi)djW - 2{2a + l){zi - Zj)w . (3.28) 

Using again Eqs. f l3.13p and f l3.20p we obtain 

+ 2wiA[{zl - zl)\s)] - 2{2a + l)«;A[(^i - Z2)\s)] . (3.29) 

Since A[(;zf - 4)|s)] = A[zf(l - 5i2)|s)], the RHS of the latter equation 
clearly belongs to Tg^-^. This shows that To(ro") C To"-^ C T^". 

Consider next the action of the operator Ti on a state of the form g{(Ji, rjv)$*^^^ , 
with g G Si~^ . From the identity 

{zidi-zi+idi+i)zi = l-^^-^^^^{l-Ki^i+i)zi + ]-{zi-zi+i){6i^i + 6i^i-i), Mi, I 
^ zi — Zi+i Z 

we easily obtain 

Ti{gzi) = {Tig)zi + 2zidig + 2ag , 
and therefore, by Eqs. (ISTIB . (ITT3|) and (IXTB|) . 

Ti =Y,Ti{gZi)\s,) = {T,g + 2g{)^^^^ + 2{ag + r^(7^)<l>(°) G T^' . 

(3.30) 

Thus TiiT-^) C Ti""^ C Ti", as claimed. 

Consider, finally, the operator T2. If g{T^_i,T]^) G ^2"^, the identity 

{zfdi-zf_^^di+i)zi = -^^-^^^{1 - Kij+i)zi + Zi{zi - zi+i){6i^i + 6i^i_i), Mi, I 

Zl — Zi^i 
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yields 

T2{gzi) = (T2g)zi + 2zfdig + Aazig , 
and hence, by Eq. fl3.17p . 

T,{g^^'^) = Y,T2{gz;)\s,) 

i 

= [T2g + 2{TN-igN-i + r^gN + 2ag)] - 2TNgN-i<^^''^ (3.31) 

clearly belongs to Tp on account of Eqs. f lS.llcp and (13.181) . The last type of 
spin states we need to study are of the form q^^~^\ where q = rArg(cri, ta?) 
with qii = 0. Since 



{zfdi - zf^^di+i)Zi ^ = ^'^'^^ (1 - K^+i)zi ^ V /, 

Zl — Zi^i 

we obtain 

T^iqzr^) = {T2q)zr^ - 2d,q + 2qzr^ , 
and thus, by Eqs. (IHTTTll and ([320]), 

r2(g$(-i)) = ^T2(g^ri)|s,) = (Tag - 2r^g^ + 2g)$(-i) - 2gi$W . (3.32) 

i 

From Eqs. (I3.11cp and (13.180 . it follows that the RHS of the previous equation 
belongs to 7^". Hence T2(7^") C 7^", which concludes the proof. □ 

Remark 3.7. We have chosen to allow a certain overlap between the different 
types of states spanning the spaces 7^". For instance, if \s) is symmetric a 
state (7(0"!, (J2, o"3)A(2;i|s)) G is also of the form /(cxi, (72, ct3)A|s). Less 
trivially, in the case of spin 1/2 the identity 

$(=^)(|s)) = ^(cri$(2)(|s)) - a2$(^)(|s))) , S,2\s) = -|s) , 

where ^'•'^^ and are respectively defined in Eqs. (13. 220 and (13.270 . implies 
that the states of the form w{ai, (T2)$*-^-* in the space T^'' can be expressed in 
terms of the other generators of this space. 

Corollary 3.8. For each e = 0,1,2, the gauge Hamiltonian leaves in- 
variant the space defined by 

'Xjn 'T"ra 'jy™ 1"^ I 1-1'^ '7~™ 
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Proof. We shall begin by showing that each space TJ^ is invariant under the 
operator J°. Note first that 

jO^O) ^ ^ jo5(2) ^ 2$(2) ^ jO J{3) ^ 3^(3) . j G Z , (3.33) 

where the states $(2) ^nd $(3) ^re defined in Eqs. flX^ . flX^ and flHITTj) . 
respectively. Using Eqs. fl3.13p . (13.151) and (13.171) one can immediately estab- 
lish the identities 

J'f = ^ifi + 2^2/2 + 3^3/3 , V/(ai, (72, ^3) , (3.34a) 

= ^1/1 + 2(T2/2 + iVr^/iv , V/((Ti, a2, r^) , (3.34b) 

J'f = ^1/1 + (iV - l)rN-ifN-i + Nr^f^ , V/(o-i, r^_i, r^) . (3.34c) 

From Eqs. (I3.33l) - p.34l) and the fact that J° is a derivation it follows that 
leaves invariant the spaces TJ^ for all e = 0,1,2. This implies that H^ 
preserves 7^" for e = 0,2, since the coefficient c_ vanishes in these cases 
(cf. Table EH]). On the other hand, for e = 1 the coefficient c_ is nonzero, 
and thus we have to consider the action of the operator J~ on the space T^. 
We now have 

j-^U) = j^U-i) ^ jGZ, (3.35) 

and, from Eq. fl3.15p . 

J-f = iV/i + 2^1/2 + rjv-i/jv , y fieri, a^, r^) . (3.36) 

Hence J~ leaves invariant the subspace of 7^" defined by the restric- 
tions fN = Qn = 0. From the obvious identity Ti(/$(°)) = (Ti/)$(°) and 
Eq. (ESQ]), together with (IXTTbl) . ([312]) and (IHIHD . it follows that the oper- 



ator Ti also preserves H". Likewise, Eqs. (I3.33P and (]3.34bp imply that 
is invariant under J°, and hence under the gauge Hamiltonian Hi. □ 

Theorem 13.61 characterizes the invariant space in terms of the subspace 
S' C S in Eq. (I3.2ip that we shall now study in detail. In fact, from the 
definition of the invariant space 7^" it follows that we can consider without 
loss of generality the quotient space where |s) ~ |s) if A(ziZ2|s)) = 

A{ziZ2\s)). For instance, from Eq. (13.190 it immediately follows that if |s) G 
S' and TT G S]^ is a permutation such that 7r(i) G {1,2} for i = 1,2, then 
7r|s) belongs to S' and is equivalent to |s). 

In the rest of this section, we shall denote |s^) simply as \s.ij) for the sake 
of conciseness. From Eq. (I3.19P it easily follows that any symmetric state 
belongs to S', since 

5^|s.,.+i) = ^^^|s), for all |s) G A(S) . (3.37) 
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On the other hand, if |s) G A(S) the corresponding state h{ai, a2)A{ziZ2\s)) 
is a trivial (factorized) state. We shall next show that the reciprocal of this 
statement is also true, up to equivalence. 

Lemma 3.9. For every |s) G S, A{ziZ2\s)) is a factorized state if and only 
if |s) ~ A|s) . 



Proof. Suppose that 



A{ziZ2\s)) = |s) ^ 



Cij Zi Zj 



is a factorized state. Since the LHS of the previous formula is symmetric, 
Cij = c for all i,j and |s) G A(S). By absorbing the constant c into |s) we 
can take c = 1 without loss of generality, and therefore 

A{ziZ2\s)) = ^ZiZj\Sij) = T2\s) =^ \Sij) = \s), i,j = l,...,N. 

i<j 

From Eq. (13.201) with f{zi, Z2) = 1 it then follows that 

A|s) = 5^|s,,) = liV(Ar-l)|s). 

i<j 

Setting I So) = |s) — A|s) and using the previous identity we obtain 

2T2 2T2 

A{ziZ2\so)) = A(2;i2;2|s)) - _ ^1^) = |s)^2 - _ -^^ A|s) = 0. 
Hence |s) ~ A|s), as claimed. □ 



By the previous lemma, it suffices to characterize the nonsymmetric states 
in S'. To this end, let us introduce the linear operator A : S — S by 

A\s) = J2Ki+i)- (3-38) 

i 

Given an element |s) = |si . . . sat) of the canonical basis of S, we shall also 
denote by {s^, . . . , s"} the set of distinct components of s = (si, . . . , sn), 
and by z/j the number of times that s* appears among the components of s. 
For instance, if |s) = |— 2, 0, 1, — 2, 1), then we can take = —2, = 0, 
= 1, so that z/i = z/3 = 2, 1^2 = 1. Consider the spin states |Xj(s)) = \xi), 
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i = 1, . . . , n, given by 



\Xi) = ^M- ^)\s"s\ . .) - Uj{uk~ 5jk)\sh^ . . .) 



(3.39a) 



\Xi 



)= J2 ^,(Nv...) + |.v...)) 



Ui = 1. 



(3.39b) 



Here we have adopted the following convention: an ellipsis inside a ket stands 
for an arbitrary ordering of the components in s not indicated explicitly. 
Note that the states (13.391) are defined only up to equivalence, and that 
IXi(s)) = Ixi(vrs)) for any permutation vr G Sat. 

Proposition 3.10. Given a basic spin state \s), the associated spin states 



Proof. Consider first a state \xi) of the type (I3.39al) . Using the definition of 
the operator A in Eq. (I3.38P we obtain 



I TrGS'iv — 2 I 

-2E E E ^.K-M^I---f^'---), (3.40) 



where the permutations tt act only on the N — 2 spin components specified 
by an ellipsis. On the other hand, we have 



N ■ N\A\s) = J2 E M^i-'^)T^\---s"s\..) 

I 7reS'jv_2 I 

+2E E E -.-i...f ...) 

I l^^mreSN-i I 

-EE E r^A^k-S,,)n\...s^s'...). (3.41) 

I l^jjfc^n 7re5Ar-2 I 
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Comparing Eqs. (13 .40^ and (13.411) we obtain 



2 (at -2 z/jjA|s) = 2(2z/, -iV)A|s). (3.42) 



This shows that any state of the form (I3.39al) belongs to S'/ ~. Suppose next 
that Ui = 1, so that \xi) is given by Eq. fl3.39bp . Since 



^I^^) = m5ZE E ^,vr(|...f.^..) + |...|V...))=4A|s), 

(3.43) 

it follows that in this case \xi) is also in S'/~ . □ 



Remark 3.11. Just as symmetric spin states, cf. Eq. (13.371) . the states \xi) 
satisfy the relation 

A\X^) = J;^MX^)■ (3.44) 
Indeed, if z/j > 1, from Eqs. (I3.39ap and (I3.42p we have 



A|X.) 



A|s) 



[uiiui -1)+N -i^,-{N - z/i)2] A|s) 
(iV_i)(2z.,-Ar)A|s) = ^^A|x.). 



On the other hand, if i/j = 1 Eqs. (13.39bp and (I3.43P imply that 

A^- 1 



A\X^) = 2( J2 '^,)A|s) = 2(Ar-l)A|s) 



Example 3.12. We shall now present all the states of the form (13.391) for 
spin 1/2. In this case, up to a permutation the basic state |s) is given by 



u N-u 



|s) = l + 



(3.45) 
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If u is either or A^, then n = 1 and thus is of the type f l3.39ap and 
proportional to |s). If z/ = 1, then n = 2 and we can take (dropping inessential 
factors) 

IXi) = ^ (|+- ■ ■ • ) + I-+ • ■ ■ )) ~ |+- • ■ ■ ) , 1X2) = I ) . 

Although the states and IX2) are linearly independent, the combination 
2|xi) + (A^ — 2) 1x2) is equivalent to a symmetric state. In the case v = N — 1 
the states \xi) are obtained from the previous ones by flipping the spins. 
Finally, ii 2 ^ v ^ N — 2 then n = 2 and the states \xi) are now given by 

IXi) = -1X2) = u{u~l)\++---)-{N-iy){N-u-l)\ ) . (3.46) 

According to the previous example, for spin 1/2 there are exactly n — 1 
independent states of the form fl3.39p associated with each basic state |s), 
up to symmetric states. We shall see next that this fact actually holds for 
arbitrary spin: 

Proposition 3.13. Let |s) be a basic spin state. Ifn is the number of distinct 
components of s, there are exactly n — 1 independent states of the form fl3.39p 
modulo symmetric states. 

Proof. Let p be the number of distinct components of s such that Vi> 1. 
A straightforward computation shows that the combination 

n n 

E 1^^) ~ (2 -P) E ^''^ - ^^^)\'''' {2-p)N{N - l)A|s) (3.47) 

is equivalent to a symmetric state. Suppose first that p ^ 2. It is immediate 
to check that in this case the set {\x-i) : « = 1, . . . , n} is linearly independent. 
If a linear combination Yll=i^i\Xi) is equivalent to a symmetric state |s), 
then |s) must be proportional to A|s), so that we can write 

n 

5^Q|x.)~A(2-p)Ar(Ar-l)A|s). 

i=l 

Hence Yl^=ii^i ~ ^)\Xi) ~ 0, and the linear independence of the states \xi) 
implies that q = A for all i. On the other hand, if p = 2 the set {\x-i) '■ 
i = 1, . . . ,n} is linearly dependent on account of Eq. f l3.47p . but removing 
one of the two states with i/j > 1 clearly yields a linearly independent set. 
It is also obvious from the coefficients of the states . . . ) that no linear 
combination XlILi '^i\Xi) can be equivalent to a nonzero symmetric state. □ 
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The next natural question to be addressed is whether the states of the 
form f l3.39p span the space E'/~ up to symmetric states: 

Proposition 3.14. The space (S'/A(S))/~ is spanned by states of the form 
(13391) . 

Proof. For conciseness, we present the proof of this result only for the case 
M = 1/2. Let S,^ denote the subspace of S spanned by basic spin states with 
z/ "+" spins, and set T,'^ = S' fl S,^. Since the operators A and A involved 
in the definition (13.211) of S' clearly preserve H^, it suffices to show that the 
states Ixj(s)) with s given by f l3.45p span the space up to symmetric 

states. Note first that the statement is trivial for u = 0,1, N — 1,N, since 
in this case the states of the form fl3.39l) obviously generate the whole space 
T,^/ ~ . Suppose, therefore, that 2 ^ u ^ N — 2, so that 

s,/~ = (|++...), |+----), I )>. 

Since the state fl3.46l) and the symmetric state (up to equivalence) 

z/(z/ - 1)1++ ■ ■ • ) + 2u{N - z/)|+ ) + (A^ - i^){N - z/ - 1)1 ) 

both belong to we need only show that (for instance) |H — ■ ■ ■ ) is not 

in i.e., that A\-\ — ■■■) is not symmetric. But this is certainly the 

case, since a state of the form 




A; = 1, 2, ... , min(z/, N — u) , 

appears in A\-\ — ■ ■ ■ ) with coefficient 2fc(z/ — 1)!(A^ — z/ — 1)! depending on 
k. □ 



3.4 The algebraic eigenfunctions 

In the previous section we have provided a detailed description of the spaces 
Ti^ C A(C[z] ® S) invariant under the corresponding gauge Hamiltonian 
H^. In this section we shall explicitly compute all the eigenfunctions of the 
restrictions of the operators to their invariant spaces This yields 
several infinit^ families of eigenfunctions for each of the models (13.21) . which 
is the main result of this paper. We shall use the term algebraic to refer 



^For the model l|3.2cp we shall see below that the number of eigenfunctions with a given 
total momentum is finite. 
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to these eigenfunctions and their corresponding energies. It is important 
to observe that the eigenfunctions of the gauge Hamiltonian if^ that can 
be constructed in this way are necessarily invariant under the whole sym- 
metric group, in spite of the fact that H^ is symmetric only under cyclic 
permutations. In fact, the explicit solutions of all known CS models with 
near-neighbors interactions (both in the scalar and spin cases) can be factor- 
ized as the product of a simple gauge factor analogous to /i times a completely 
symmetric function [10,48,75,124]. This, however, does not rule out the ex- 
istence of other eigenfunctions of the gauge Hamiltonian if^ invariant only 
under the subgroup of cyclic permutations, which is indeed an interesting 
open problem. 

Case a 

We shall begin with the model fl3.2al) . which is probably the most interesting 
one due to the rich structure of its associated invariant flag. In order to find 
the algebraic energies of the model, note first that one can clearly construct 
a basis Bq of TYq whose elements are homogeneous polynomials in z with 
coefficients in S. If / G Bq has degree k, then J^f = kf and To/ has degree 
at most k — 2. li Bq is ordered according to the degree, the operator Hq 
is represented in this basis by a triangular matrix with diagonal elements 
Eq + kco, where k = 0, . . . ,n is the degree. Thus the algebraic eigenfunctions 
are the numbers 

Ek = Eo + 2ku, A; = 0, 1, ... . 

We shall next show that the algebraic eigenfunctions of Hq can be ex- 
pressed in closed form in terms of generalized Laguerre and Jacobi polyno- 
mials. The computation is basically a two-step procedure. In the first place, 
one encodes the eigenvalue problem in the invariant space Tig as a system of 
linear partial differential equations. The second step then consists in finding 
the polynomial solutions of this system. 

Regarding the first step, we shall need the following preliminary lemma: 
Lemma 3.15. The operator Hq preserves the following subspaces o/T-Cq: 

(3.48) 

^i2|s) = |s), (3.49) 
5i2|s) = -|s), (3.50) 

where f, g, h, h, w are as in the definition ofT^ in Theorem \3.6\ and ^^^^ 
$(2)^ $(3) respectively given by (IX^ . ( KT^i and (IXTTj) . 



^o,|s)=^o,|s) + P^'^), 
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Proof. The identity To(/<l>(°)) = (To/)<l>(°) and Eqs. (ESI, ^131241) and (Km 
clearly imply that the subspace '^^o is) is invariant under Hq. Consider next 

the action of Hq on a function of the form h^^'^\ Since |s) is symmetric 
under 512, we can replace A[(zi + -22) |s)] by 2$^^^ in Eq. (I3.26p . Secondly, 
any state |s) G S' satisfies the identity 

El<m) = ]^A|s). (3.51) 

i 

Indeed, we already know that the previous identity holds for symmetric states 
(cf. Eq. (E37D) and for states of the form fICTD (cf. Eqs. (1^^^ and flCTD l 
On the other hand, by Proposition 13.141 every state in E' is a linear combi- 
nation of a symmetric state, states of the form (13.391) . and a state |s) such 
that A(ziZ2|s)) =0. But for the latter "null" state |sjj) = for all i < j, 
and hence A\s) = A|s) = 0. Therefore, Eq. (I3.26P can be written as 

To = ^Toh + Sh)^^^^ + 4/ii$(i) - . (3.52) 

From the previous equation and Eq. (lOHll it follows that Ho{h(^^'^^) e H^.y 
Finally, if 5i2|s) = -|s), Eq. fl3:29ll reduces to 

To(w$(3)) = (Tqw + 12^2)$^^^ + 4wi<l>(2) - 4(2a + 1)^$^^) , (3.53) 

which, together with Eq. (l333ll . implies that Ho{w^^^^) G H^^^y □ 

Remark 3.16. The requirement that |s) be symmetric (respectively antisym- 
metric) under 5*12 in the definition of the space '^^o [s) (respectively |g^) is 
no real restriction, since the antisymmetric (respectively symmetric) part of 
|s) does not contribute to the state <l>*^^^ (respectively ^^^^). 

By the previous lemma, we can consider without loss of generality eigen- 
functions of Hq of the form 

$ = + ^$(1) + /i<|.(2) + /^${2) + y^${3) ^ deg$ = fc, (3.54) 

where the spin functions $(2) and are all built from the same spin 
state |s). Note that we can assume that hw = 0, and that the spin state |s) is 
symmetric under S12 and belongs to S' if /i 7^ 0, whereas it is antisymmetric 
under Su if w 7^ 0. 

Using Eqs. (MD, (IXg^ and it is straightforward 

to show that the eigenvalue equation Hq^ = {Eq + 2ku!)^ is equivalent to 
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the system 

[-To + 2cj( J° + 3-k)]w- I2w2 = , (3.55a) 

[-To + 2uj{J° + 2 - k)]h - 8h2 = , (3.55b) 

[-To + 2cj( J° + 2-k)]h-8h2 = 6^3 + Awi , (3.55c) 

[-To + 2cj( J° + l-k)]g-4:g2 = 4(/ii + h) - 4(2a + l)w , (3.55d) 

[-To + 2^( jO -k)]f = 2 ((71 + (2a + l)h - h) . (3.55e) 

Since / and g are linear in CT3, we can write 

/ = P + o"3g, g = u + a3V, (3.56) 



where p, g, m and v are polynomials in cti and CT2. Taking into account 
that the action of Tq on scalar symmetric functions is given by the RHS of 
Eq. (13121) with e = 0, and using Eqs. (EHi), (Km and (KM^ . we finally 
obtain the following linear system of PDEs: 



[Lo 


~2uj{k- 


3) 


~\w - 


- \2W2 


= 0, 


(3.57a) 


[Lo 


- 2uo{k- 


2) 


]h- 


8h2 = 


0, 


(3.57b) 


[Lo 


-2uj{k- 


2) 


]h- 


8h2 = 


6v + iwi , 


(3.57c) 


[Lo 


-2u{k- 


1) 


]u- 


4M2 = 


4/11 + Ahi + 6cr2f 1 














+ 6(2a + l)(Tif-4(2a + l)w, 


(3.57d) 


[Lo 


- 2uj{k- 


4) 


]v- 


16V2 = 


= 0, 


(3.57e) 



(To - 2ujk)p = 2ui + 2(2a + l)h - — — - h + Gaagi + 6(2a + l)aiq , 

(3.57f) 

[Lo -2u;{k- 3)] q - 12^2 = 2vi , (3.57g) 
where 

To = -{Ndl^ + Aaid^.d^, + 4a2dl^ + 2{2a + l)Nd^,) 

+ 2uj{aid^, + 2a2d^,) . (3.58) 

As a consequence of the general discussion of the previous Section, the lat- 
ter system is guaranteed to possess polynomials solutions. In fact, these 
polynomial solutions can be expressed in closed form in terms of generalized 
Laguerre polynomials and Jacobi polynomials pIP'^\ 
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Theorem 3.17. Let 



2r^ 



-1 



a^N(a+^^-^, /3 = /3(m) = l-m-7v(a+^), t 

where x — Xi is the center of mass coordinate. The Hamiltonian 

Hq possesses the following families of spin eigenfunctions with eigenvalue 
Ei^ — Eq + 2uj{2l + m), with I ^ and m as indicated in each case: 



(0) _ 



= ^x-^-^Ll\ujr')P^^^\t) ($(1) - X $W) , m ^ 1 , 



Im 



m ^ 2, 



(2) 



(^Plmj_i W + ^2Lir^j + _ 1) J-i 



m ^ 2, 
m ^ 3 , 



PL^^f(t)($(')-2a;$(2)) 



+ 2X I P|^ m,-3 J 



2[fJ-l m 



2r 



J- p(a+3,/3) 

3 L^J 



^ p(«+2,/3) 

2[fJ -1 L^J 



(^) 



+ (1 - (-i)"^) o^^r^^^^ -pSf(^) )$^°^ 



*S = l^x-%'\ 



Lor 



+ 



2m(m-2) L^J 

,(a+3,/3) 
LfJ-2 



m ^ 3 , 



2 
/ 1 



L2(L^J , „ 

3 



m ^ 4 . 



Here 



= A(a;^|s)), = A(xia;2|s)), S^^^ = A{xiX2{xi - X2)\s)), 
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where the spin state |s) is symmetric under Su and belongs to S' for the 

y(2) 

" Im ' 



eigenfunction and is antisymmetric under Su for the eigenfunction 



^ (s) 

The functions ifm, 4>m Ci^^d Xm O'fe polynomials given explicitly by 



m + 2a + 2 1^+2,(3-2) _ p(a+3,/3-l) _ 4a + 7 („+2,/3-l) 1 p(a+3,/3) 

1 rn -L rn -1 -L rn -1 ~r „ ± rri_ o 5 

777,-1 2 2^ 777_l2^ 32^ 

, _ p(a+4,/3-l) _ ^ p(a+3,/3-l) _ m + 2q; + 3 (a+2,/3-1) 

2 ^ 2 ^ (m — l)(m — 3) 2 ^ 

_ 1 p{a+4,/3) m + 2a - 1 (q+3,/3) 

3 2^ m — 3 2 ^ 
_ 3m + 2a (a+2,/3-1) 2m - 7 /a+3,i3-i) _ (a+4,/3-1) 
(m— lj(m — 3) 2 ^ m — 3 2 ^ 2 ^ 

m + 2a + 2 (a+2,/3) _ m + 2a (a+3,i3) 1 p(a+4,/3) 



(m — l)(m — 3) ^~ m — 3 3 

/or even m, and 

,„ _ 9p(Q+2,/3-l) _ p(a+3„;3-l) 1 p(a+3,/3) m + 2a + 2 

V-^m — m-l m-1 T J m-3 "T / ^\ m-3 

3 m(m — 2) 

_ m + 2a + 2 («+2,/3) 

2 m — 3 5 

, _ p(a+4,/3-l) _ 2m - 5 (a+3,/3) _ 1 p(a+4,/3) m + 2a - 1 (a+3,/3) 

'rm m — 3 „ m — 3 „ m-5 "T m — 5 ; 

m — 2 3 m — 2 

_ 2m - 3 („+2,/3-i) 2(m - 3) (a+3,/3-1) _ (a+4,/3-i) 

Xm I „\ ^ m-3 ~r 1 m-3 

m(m — 2) m — 2 

_ m + 2a + 1 (a+2,/3) _ r» + 2a (a+3,/3) 1 p{a+4,/3) 

m(m — 2) m — 2 3 

/or odd m. For every n = 0,1, . . . , the above eigenf unctions with 21 + m ^ n 
span the whole HQ-invariant space fiH^. 



Proof. Recall, to begin with, that the algebraic eigenfunctions of Hq are of the 
form \E' = with fi given in Table [3TT] and $ an eigenfunction of Hq of the 
form fl3.54l) - fl3.56p . In order to determine $, we must find the most general 
polynomial solution of the inhomogeneous linear system fl3.57l) . From the 
structure of this system it follows that there are seven types of independent 
solutions, characterized by the vanishing of certain subsets of the unknown 
functions p, q, u, v, h, h, w. These types are listed in Table 13.21 where in the 
last column we have indicated the eigenfunction of Hq obtained from each 
case. We shall present here in detail the solution of the system f l3.57p for the 
case q = v = h = h = w = and m 7^ 0, which yields the eigenfunctions of 
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Table 3.2: The seven types of solutions of the system (13.571) and their corre- 
sponding eigenfunctions. 









Conditions 


Eigenf unct ion 


Q = 


u = 


V = 


h = h = w = 0, PT^O 


T (0) 

Im 


u = 


V = 


h = 


h = w = 0, q 




q = 


V = 


h = 


h = w = 0, u 


^ Im 


Q = 


V = 


h = 


w = 0, hy^O 




Q = 


V = 


h = 


w = 0, h^O 




Q = 


V = 


h = 


0, Wy^O 




h = 


w = 


= 0, 







the form \E'jJ (the procedure for the other cases is essentially the same). In 
this case the system (I3.57P reduces to 

[Lq -2uj{k- 1)] M - 4u2 = , (Lo - 2ujk)p = 2ui . (3.59) 

Let us begin with the homogeneous equation for u. We shall look for poly- 
nomial solutions of this equation of the form u = Q{o'i,cr2)R{o'2), where Q 
is a homogeneous poljTiomial of degree m — 1 in z and i? is a polynomial of 
degree / in cr2, so that k = deg$ = 2/ + m by Eq. (I3.54p . From Eq. (I3.58P 
and the homogeneity of Q we have 

LoiQR) = {LoQ)R + Q{LqR) - 4aiQiR2 - 8a2Q2i?2 
= {LoQ)R + Q{Lo - 4(m - l)d^,)R . 

Hence the equation for u can be written as 

R{Lo - 49,,)g = Q{-Lo + Amd^, + 41uj)R, 

where Lq = Lo|a;=o- Since (Lq — 4(9o-2)Q is a homogeneous polynomial of 
degree m — 3 in z, both sides of the latter equation must vanish separately. 
We are thus led to the following decoupled equations for Q and R: 

{Lo-Ad^,)Q = 0, (3.60) 
{- Lo + 4:md^^+Alu)R = 0. (3.61) 

In terms of the variable p = uja2, Eq. (I3.6ip can be written as 

4cj£-^(/2) = , 
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where 

Ct = pdl + + p)dp + v (3.62) 

is the generahzed Laguerre operator. Hence R is proportional to the gener- 
ahzed Laguerre polynomial LJ^{uja2)- On the other hand, we can write Q = 
a^^^P{t) where P is a polynomial in the homogeneous variable t = ^^2. _ x 
With this substitution, Eq. fl3.60p becomes 



L 2 



r(7,<5) 



(P) = 0, 



where the Jacobi operator J'u is given by 

j-M) = (1 _ ^2)^2 ^ _ ^ _ ^ ^ ^ 2)t] dt + iy{iy + ^ + 5 + l). 



Thus P{t) is proportional to the Jacobi polynomial Pj' Ztl fHt) , so that we 

L 2 J 

can take 



u 



a 



m—l 
1 



L7^{uja,)P, 



it). 



(3.63) 



We must next find a particular solution of the inhomogeneous equation 
for p in 03.591) . since the general solution of the corresponding homogeneous 
equation yields an eigenfunction of the simpler type \E'[|)]. Since 



Ui 



[ra 



l)PL^^f (t)-2(t + l)P^(^f(t) 



(where the dot denotes derivative with respect to t), we make the ansatz 
p = Q{(Ti,a2)R{cr2), where Q is a homogeneous polynomial of degree m in 
z and P is a polynomial of degree / in a2- Substituting this ansatz into the 
second equation in (13.591) and proceeding as before we immediately obtain 

R(LqQ) +Q{Lo- Amd^^ - 41uj)R = 2ui . 

If we set R = L^'^{uja2) the second term of the LHS vanishes, and canceling 
the common factor LJ^{uj(T2) we are left with the following equation for Q: 



LqQ = 2a" 



m-2 



[m - l)p(^f (t) - 2(t + l)P,^^'|^^(t) 



p(a+l,/3), 



The form of the RHS of this equation suggests the ansatz Q = a^P(t), with 
P a polynomial in the variable t. The previous equation then yields 



1 

2iV 



[m-l)Pf^f (t)-2(t 



DP 



it) 



(3.64) 
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From the definition of tlie Jacobi operator we easily obtain 



(«,/3) 

LfJ 



("+l,/3) 



L^wj +il + t)dt- -{m-l] 



wliicli implies that P = — — P^Ztlf^ is a particular solution of Eq. (13.641) . 

N L 2 J 

Hence 



1 



(3.65) 



is a particular solution of the inhomogeneous equation in (I3.59p . We have 
thus shown that $ = p<|)W + u^^^\ with u and p respectively given by 
Eqs. (I3.63P and (I3.65p . is an eigenfunction of Hq with eigenvalue Eq + 2u{21 + 
m). Multiplying $ by the gauge factor /x we obtain the eigenfunction of 
Hq in the statement. 

It remains to show that the states (fc = 0, . . . , 4), \&[^ and \E';^ with 
2/ + m ^ n generate the spaces (I3.48p - p.50p . Consider first the "monomials" 
of the form /iO"5^(T2 which belong to /uH^,^^ if 21 + m ^ n. We can 

order such monomials as follows: we say that 

21' + m' < 2l + m, or 21' + m' = 2l + m and m' < m. From the expansion [96, 
Eq. 8.962.1] 

Pl'''\t) = ^ E ^ (-^)^(7 + S + U + 1).(7 + k + 1)._, (1 - , 

where {x)k is the Pochhammer symbol 

{x)k = x{x + 1) ■ ■ ■ (x + A; - 1) , 

it follows that p!)'^'^\o) > provided that 7 + 1 > and 7 + 5 + 2z/ < 0. In 
particular, Pf^f (0) > smce 



1 



1 



a+l = N{a+-]--> N-->0, a+p+2 



m 




m 


1 




= 2 




— m — 


y 




_ 2"_ 


2 



1 

2 



Hence we can write 



<)=/x<l>W(Q^ar4 + l.o.t.) 



where 7^ 0, so that 

: 2/ + m ^ n) = (/lalX"^^"^ : 2/ + m ^ n) . 
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Likewise, a similar argument shows that for m ^ 1 

Ux"'-^L;'^{iur^)P^^[^\t) <l>(^) : 2l+m ^ n) = (/i<-V^$(^) : 2l+m ^ n) , 

L 2 J 

and therefore 

: 2/ + m < n) = /i<-V^$(i) : 2/ + m ^ n) . 

Proceeding in the same way with the remaining spin eigenfunctions we can 
finally show that 

(vi>;^^ = 0, . . . , 4 , 2/ + m ^ n> = , 

and that 

/^^o,|s) + : 2/ + m ^ n> = /.TY^^I,^ , |s) G S', ^i^ls) = |s), 

/^^o,|s) + : 2/ + m ^ n> = /.T^^,,,^ , S^^ls) = -|s), 

as claimed. □ 
Case b 

Since Cq = 4u; in this case, reasoning as before we conclude that the algebraic 
energies are the numbers 

Ek = Eq + Akuo, A; = 0, 1, ... , 

where k is the degree in z of the corresponding eigenfunctions of if i. We shall 
see below that these eigenfunctions can be written in terms of generalized 
Laguerre polynomials. To this end, we begin by identifying certain subspaces 
of invariant under Hi. 

Lemma 3.18. For any given spin state |s) G S, the operator Hi preserves 
the subspace 

K\s) = (/(^i' ^2) gicri) $(1) : f22 = 0) C (3.66) 

where f and g are polynomials of degrees at most n and n—1 in z, respectively, 
and is given by ( I3.22p . 

Proof. The statement follows from the obvious identity Ti (/$''°'') = (Ti/)$*^'') 
and Eqs. (imb|) . (IXT^ . (IHTTUD . flCTj) . (IX^ and (1S3SD. □ 
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By the previous lemma we can assume that the eigenfunctions of Hi in 
are of the form 



)(!) 



deg $ = A; ^ n . 



(3.67) 



From Eqs. fl3.30p . 03.331) and fl3.35p it easily follows that the eigenvalue equa- 
tion = {Eq + Akuj)^ can be cast into the system 



1 



9 - 2^1 



-Ti +cj(J°+ 1 - k) 

-Ti+uj{J^ -k)- {^ + l^J']f = i^ci + h 



0, 



Since / is linear in a2 (cf. Eq. f l3.66p ). we can write 

f = p + (Tag, 



(3.68a) 
(3.68b) 

(3.69) 



where p and q are polynomials in ai. Using Eqs. (13.11bp . fl3.12p . fl3.16p . 
(I3.34bp and (13.360 we can easily rewrite the system (13.680 as follows: 



[Li-ujik-l)]g-2gi = 0, 
[Li - uj{k - 2)]q - Aqi = , 

{Li - ujk)p = (2a + b + ^^g + 2(Aa + b + ^^aiq . 



where 



Li = -aidl^ + 



uai- {2a + b + ^]N 



da 



(3.70a) 
(3.70b) 

(3.70c) 
(3.71) 



The last step is to construct the polynomials solutions of the system fl3.70p . 
which can be expressed in terms of generalized Laguerre polynomials, ac- 
cording to the following theorem. 

Theorem 3.19. The Hamiltonian Hi possesses the following families of spin 
eigenfunctions with eigenvalue Ek = Eq + Aku: 



(0) 



a-l, 



$(0) 



k >0 



^« = f,L"^+l (ur^) [iV$(i) - r2$(0)] , k^l 



(2) 



a+3| 



Ar(a + l)^xf-/?r^j k^2, 



wherea = N{2a+b+l), (3 = A^(4a+6+|), and$(j) = A(xi>)), with j = 0,1 
and |s) G S. For each |s) G S and n = 0, 1, ... , the above eigenfunctions 
with k span the whole Hi-invariant space fiTC^^^y 
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Proof. As in the previous case, the algebraic eigenfunctions of Hi are of the 
form \l/ = where /i is given in Table 13.11 and $ is an eigenfunction of 
Hi of the form fl3.67l) - fl3.69l) . The functions p, q and g are polynomials in ai 
determined by the system (13.701) . which in terms of the variable t = ujai = 
ujr"^ can be written as 



i^T-l9 = = , Cr'p =-^9-^,tq, (3.72) 



where is the generalized Laguerre operator (cf. Eq. (13.621) ). The general 
polynomial solutions of the first two equations in (13.721) are respectively given 
by 

g = CiLtll{t), q = c^Lltl{t). (3.73) 
On the other hand, from the elementary identity 



Cl{t^Lltf{t))=l{l + \)t'-'Lltf{t) 



it follows that the general polynomial solution of the third equation in (I3.72p 
is given by 

Equations (I3.73P and (13.741) immediately yield the formulas of the eigenfunc- 
tions in this case. The last assertion in the statement of the theorem follows 
from the fact that the functions p, q and g in Eqs. (13.731) and (13.741) are the 
most general polynomial solution of the system (13.721) . □ 



Remark 3.20. It should be noted that for = the potentials (I3.2al) and fl3.2bp 



scale as under dilations of the coordinates (as is the case for the orig- 
inal Calogero model). The standard argument used in the solution of the 
Calogero model shows that there is a basis of eigenfunctions of these models 
of the form /i(x)L^(c<jr^)F(x), where F is a homogeneous spin- valued func- 
tion. The eigenfunctions presented in Theorems 13.171 and 13.191 are indeed of 
this form. 



Case c 

In this case the number of independent algebraic eigenfunctions is essentially 
finite. We shall take, for definiteness, the plus sign in the change of variable 
listed in Table 13.11 (it will be apparent from the discussion that follows that 
the minus sign does not yield additional solutions). 
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Let us first note that the potential for this model is translationally invari- 
ant, so that the total momentum operator P = —i^f.^x^. commutes with 
the Hamiltonian H2. Hence the eigenfunctions of H2 can be assumed to have 
well-defined total momentum. Equivalently, since 

= 2J°, (3.75) 

^k = -^ log 2fe 

the eigenfunctions of H2 can be assumed to be homogeneous in z. Let $ be a 
homogeneous eigenfunction of H2 of degree k and eigenvalue E, so that \E' = 
/i$ is an eigenfunction of H2 with total momentum 2k (cf. Eq. (13.751) ) and 
energy E. By Eq. (13.751) . the function t^'^ clearly has total momentum 2{k + 
NX). In fact, the following lemma implies that rj^^^ is also an eigenfunction 
of H2 with a suitably boosted energy: 

Lemma 3.21. Let ^ be a homogeneous eigenfunction of H2 of degree k 
and eigenvalue E. Then is an eigenfunction of H2 with eigenvalue 
E + 8kX + ANX^. 

Proof. From the identity 

I I 

where Di = Zidi, we immediately obtain the following expression for the 
gauge Hamiltonian H2'. 

-{H2-E,) = J2D^ + aJ2'-^±^iD,-D.,^,) 

4 Zi — Zi+i 

_2^^^^£m ). (3-76) 

Since t^^DiT'^ = Di + X for any real A, it follows that 

t^^H2Tn =H2 + 8AJ° + 4NX^ . 
Taking into account that = A;$, we conclude that 

H2 (r^$) = {E + 8kX + ANX^) (r^<l>) , 
as claimed. □ 

By the previous discussion, in what follows any two eigenfunctions of H2 
that differ by a power of shall be considered equivalent. From Theorem l3.6l 
and Corollary l3.8l it easily follows that in this case the number of independent 
algebraic eigenfunctions is finite, up to equivalence. More precisely: 
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Lemma 3.22. Up to equivalence, the algebraic eigenf unctions of H2 can he 
assumed to belong to a space of the form 

H2,|s> = (f^i, (TiTN-i, r^) $ + (1, r;v-i) + (1, ^i) r^^^'^) (3.77) 

for some spin state |s), where ^^'^'^ is given by fl3.22p . 

Proof. Given a spin state |s), the obvious identity 

T2(/$W) = (T2/)<l>(°) (3.78) 

and Eqs. (13.71) . (13.311) . (I3.32p . and (13.331) imply that the gauge Hamiltonian 
H2 preserves the space 

'Hl\.^ = {f^'^'\g^'^'\Tm^^-'^). (3.79) 

where /, g and q are as in the definition of 7^" in Theorem l3.6[ Let $ G ^2 |s) 
be an eigenfunction of H2, which as explained above can be taken as a 
homogeneous function of z. From the conditions satisfied by the functions 
/, g and q in (13.791) and the homogeneity of $, it readily follows that $ G 
'Tn'^2,\s) for some /, as claimed. □ 

Theorem 3.23. The Hamiltonian H2 possesses the following spin eigenfunc- 
tions with zero momentum 



^3 = 



2a 



2a +_ 



_<|,(o) + ^cos(2(xi -Xj))|sj) , "$4 = f^'^sm{2{xi- Xj))\sj), 



where x is the center of mass coordinate. Their energies are respectively given 
by 

Eo, Ei,2 = ^o + 4(2a + l-i^), E3,4 = ^o + 8(2a + l). 

Any algebraic eigenfunction with well-defined total momentum is equivalent 
to a linear combination of the above eigenf unctions. 

Proof. By Lemma I3.22[ in order to compute the algebraic eigenf unctions of 
H2 it suffices to diagonalize H2 in the spaces ?^2,|s> given by (13.771) . Equa- 
tions (13.311) . (13.321) . (I3.33P and (13.780 . and the fact that H2 preserves the 
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degree of homogeneity, imply that the following subspaces of 7^2, |s) are in- 
variant under H2'- 

(alr;v-l,r^)$^°^ (3.80a) 
(ai<l>(°),<l>«), (cTir;v-i$^°\r,v$^°\r^-i$^'^) , (3.80b) 
(r^_i<l>W), (r^_i$(°\rAr$(-i)), (airA._i$(°), r^$(°), air^^^'^)) . (3.80c) 

From Eqs. fl3.80p it follows that the alternative change of variables = 
g-2ixfc (^QQ^ j^Q^ yield additional eigenf unctions of H2. Indeed, the latter 
change corresponds to the mapping Zk ^ 1/zk, which up to equivalence 
leaves the subspace f l3.80ap invariant and exchanges each subspace in ( ]3.80bll 
with the corresponding one in fl3.80cp . For this reason, we can safely ignore 
the subspaces fl3.80cp in the computation that follows, provided that we add 
to the eigenfunctions of H2 obtained from the subspaces (I3.80bp their images 
under the mapping Zk ^ 1/zk- 

For the subspaces (raiaD - flSml) . using Eqs. (fSIID, dSH, flH3g|)_flH3TD . 
fl3.33p . and fl3.34cp we easily obtain the following eigenfunctions of if 2: 

tn^^°\ E = Eo + AN, (3.81a) 

E = Eo + 4(2a+ 1), (3.81b) 

TN-i'^^^^ - E = Eo + 4{N + 4a + 2). (3.81c) 

We have omitted the two additional eigenfunctions 

(^'^--|tt)*'°' 

from the above list, since they are respectively obtained from ( 13.81bp and f l3.81cp 
when the spin state |s) is symmetric. The eigenfunctions (13.810 are equivalent 
to the following "zero momentum" eigenfunctions: 

E = Eo, (3.82a) 
r-'/^$«, E = Eo + A{2a+l-j^'), (3.82b) 

^<l>«--^<l>(°), E = Eo + 8(2a+l), (3.82c) 

tn 2a + 1 

where the energies have been computed from those in Eqs. (13.810 using 
Lemma I3.21[ The eigenfunctions of H2 listed in the statement are readily 
obtained from these spin functions together with the transforms of ( 13.82bp 
and f l3.82cO under the mapping Zk ^ l/zk- □ 
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Remark 3.24. If the spin state |s) is symmetric, then |sj) = ^^^^/N for all 
i, and one easily obtains from Theorem 13.231 the following eigenfunctions of 
the scalar Hamiltonian ifn'^: 



1,2 



/i 



cos 
sin 



/i 



aN 



2a + 1 ' ^ 

Kj 



f> cos(2{xi—Xj)') 
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Chapter 4 

NN spin chains of HS type 



4.1 Introduction 

As we saw in Chapter [21 one can successfully solve spin chains of HS type by 
means of Polychronakos's freezing trick, which essentially consists in comput- 
ing the large coupling constant limit of the quotient of the partition functions 
of a spin CS model and its scalar counterpart. 

In the present chapter we shall extend this technique to the case of QES 
models of CS type with spin. This extension will prove crucial in the analysis 
of the spectrum of the novel A^-site chain 

y^ = J2(^i-^i+irHi-su+,), (4.1) 

i 

with sites ^ G Cq C defined by the algebraic system 

^^= , \ + , \ > l^t^N. (4.2) 

Here and in what follows, we identify = ^n, (,n+i = (.i (and similarly for 
Xi and Si^i+i), and keep denoting by Cq the domain fl3.3ap and by M the 
particles' spin. It is not difficult to check that the spin chain (14.11) . which 
we shall simply call the NN chain, is related to the NN Hamiltonian Ho in 
Eq. (I3.2al) along the lines of the freezing trick. 

We shall see in this chapter that this spin chain is in some sense in- 
termediate between the Heisenberg chain (short-range, position-independent 
interactions) and the PF chain (long-range, position-dependent interactions), 
from which one can derive the Hamiltonian (14.11) by retaining only nearest- 
neighbors interactions. Moreover, H turns out to be QES in the sense that a 
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nontrivial proper subset of the spectrum can be explicitly computed for any 
value of N and M. While QES models have long played a relevant role in 
the theory of CS systems [179, 180], this is to our best knowledge the first 
QES phenomenon spotted in the (as has been made clear in this dissertation) 
closely related field of HS spin chains. 

The results in this chapter can be easily modified to deal with the NN 
chains 

H2 = ^ sin"^(^i - ^i+i) (1 - Si^i+i) , 

i 

where the chain sites are respectively defined by 

^. = 1 + ^ + ^ 

S« C C2 _ f 2 £-2 _ rr2 ' 

Si Si Sj+1 Si Sj-1 

= cot(^i - ^i+i) + cot(^i - . 

These spin chains are obtained from the NN models (13.2bl) and fl3.2cp re- 
spectively by a freezing trick argument. Note, however, that the sites of the 
second chain turn out to be equally spaced, so H2 is essentially the Heisenberg 
chain. 

This chapter is organized as follows. In Section 14^2] we analyze the distri- 
bution of the chain sites using both exact results and numerical calculations. 
In Section we show the QES nature of the chain by providing several fam- 
ilies of eigenvectors, valid for any choice of N and M. The main technical 
tool is a nontrivial refinement of the freezing trick developed in Chapter [2] 
which also applies to QES systems. Finally, we discuss some qualitative 
properties of the spectrum and formulate an intriguing conjecture claiming 
that an eigenvalue of the chain is algebraic if and only if it is an integer. 

The material presented in this chapter is taken from Ref. [63]. 

4.2 The chain sites 

We shall begin with a detailed analysis of the distribution of the chain sites on 
the line in which we shall combine exact results and numerical computations 
to extract as much information as possible. An easy observation, which stems 
from the spin chain's connection with the dynamical model fl3.2al) that will 
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be developed in the next section, is that the sites are given by the coordinates 
of a critical point of the function 



-^W = J^loglxj - — , (4.3) 



which is the scaled logarithm of the ground state function of the NN Hamil- 
tonian fl3.2al) . 

The first questions to ascertain refer to the existence and uniqueness of 
the critical points of A. In this direction we have the following result, which 
is roughly analogous to Proposition (12.201) : 

Proposition 4.1. The function A has a unique critical point ^ in Cq, which 
is a hyperbolic maximum. Moreover, ^ satisfies 

6 = -e7v-m, 2 = l,...,iV. (4.4) 

Proof. The existence of a maximum of A in Co is clear, since it is continuous 
in Co and tends to — oo both on its boundary and as r ^ oo. Uniqueness and 
hyperbolicity follow from the fact that the Hessian of A is negative definite 
in Co- Indeed, by Gerschgorin's theorem [96, 15.814], the eigenvalues of the 
Hessian of A at x lie in the union of the intervals 



I oxf oxf 



where 7j = 



dxidxj 



Since 



[Xi-Xi±i) , — — — = if J 7^ ± 1, 



OXiOXi±i OXiOXj 

we have 

oxf 

and thus all the eigenvalues of the Hessian of A are strictly negative. 
Furthermore, the mapping 

Xi -XN-i+i , i = l,...,N, 

is a symmetry of A and maps the domain Co into itself. Since ^ is the unique 
critical point of A in Cq, it must be a fixed point of the above transformation, 
and thus Eqs. (lO) follow. □ 



96 



CHAPTER 4. NN SPIN CHAINS OF HS TYPE 



From the previous proposition (and also directly from Eq. (14. 21) ). it im- 
mediately follows that the center of mass of the spins vanishes, i.e., 

J2^. = 0. (4.5) 

i 

Proposition 4.2. = A^. 

Proof. Let us write A in terms of the variables r = ||x|| G and y = x/r G 
§^-1, i.e., 

A = log IVi - + A^logr - y . 

i 

Obviously 

as we wanted to show. □ 

It is convenient to think of the spins as if located on a circle. In this 
way, the chain (14. ip truly presents only near-neighbors interactions in spite 
of the fact that the first spin interacts with the last one. Indeed, if the site 
coordinate C,i is understood as an arc length in a circle of radius 

tn = , (4.6) 

71 

as shown in Fig. 14. 11 the spins at the sites .^i and become nearest-neighbors 
and, moreover, the strength of the interaction is inversely proportional to the 
squared distance between the spins, measured along the arc. 

We saw in Section [L8] that the site coordinates of the PF chain define the 
A^-th Hermite polynomial. It is natural to wonder whether a similar result 
holds for the family of monic polynomials 

PN{t) = l[{t-^,), ^gCoCM^, (4.7) 

i 

defined by H, which can be regarded as the closest NN analog of the PF 
chain. In the following proposition we answer this question negatively: 

Proposition 4.3. The polynomials (14. 7p do not form an orthogonal family. 

Proof. In order to show that the family {pn} is not orthogonal, it suffices 
to prove that the polynomials (14. 7p do not satisfy a three-term recursion 
relation of the form 



PN+l{t) = tPN{t) + CNPN^lit) , 
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with Ctv constant. In fact, solving the system fl4.2p for = 2,3,4 one 
immediately obtains 

P2{t) =t'-l, ps{t) = - f , P4{t) = - 2t2 + 1 . 

Thus the previous recursion relation already fails for = 3. □ 

The sites are certainly not equally spaced; in fact, our numerical compu- 
tations lead to the following 

Claim 4.4. For N ^ 1, the site distribution follows the Gaussian law with 
zero mean and unit variance. 

Note that, should the site distribution follow a Gaussian law, its mean 
would be zero on account of Corollary 14. 5[ We have solved numerically 
Eq. (14.21) for the positions of the chain sites for up to = 200 spins, and 
the results completely support the validity of the above conjecture. More 
precisely, we claim that the cumulative density of sites 

J^{x) = N-'Y.9{x-^i), (4.8) 

i 

where 9 is Heaviside's step function, is asymptotically given by 

F{x) = ^ 1 + erf {x/V2) . (4.9) 



98 



CHAPTER 4. NN SPIN CHAINS OF HS TYPE 



T{x),F{x) 




Figure 4.2: Cumulative distribution functions T{x) (at its dis- 
continuity points) and F{x) (continuous grey line) for = 150 
spins. 



The agreement between the functions JF and F is remarkably good for > 
100 (see Fig. 14.21 for the case = 150) and increases steadily with N . 
For example, the mean square error of the fit for 100, 150 and 200 spins 
is respectively 2.6 x 10-^ 1.1 x 10"^ and 7.9 x 10"^ 



Heuristic proof of Claim \4-4\ The fact that for large N the cumulative den- 
sity of sites is well approximated by the Gaussian law 04.91) can be justified 
by the following semi- rigorous argument. Let x(t, A^) be a smooth function 
such that x{i, N) = for i = 1, . . . ,N, and define the rescaled function 
y{s,e) = x{s/e, 1/e). By Eq. (14. 2p . the latter function must satisfy the rela- 
tion ^ ^ 

— ^ — 7 T + — ^ — 7 7 = yi^^ e) (4.10) 

y{s,e) - y{s - e,e) y{s, e) - y{s + e, e) 

for e = 1/A^ <^ 1 and s = j^, j^, . . . ,1. Let us now assume that Eq. (I4.10p 
holds for all s G M and all e <^ 1. Writing 



1/(5, ej 



k=0 



and using the expansion 

yis, e) - y{s ±e,e) = T y',{s)e - ± y[{s)) + 0{e 
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the leading term in Eq. fl4.10p yields the differential equation 

Vo = Z/o % • 

The general solution of this equation is implicitly given by 

s = co + cied{yois)/V2) . (4.11) 

Hence, up to terms of order e = 1/A^, the cumulative distribution function 
of the chain sites (normalized to unity) is approximated by the continuous 
function 

F{x) = Co + Ci erf (x/ V^) . 

The normalization conditions -F(— oo) = and F{oo) = 1 imply that Cq = 
Ci = 1/2, and thus the empiric law (14. 9 p is recovered. 

From Eq. (14. lip (with cq = Ci = 1/2) it follows that the site can be 
determined up to terms of order 1/A^ by the formula 

6 - V2 erf-^ (^^) . (4.12) 

Note that, were the sites exactly given by the previous formula, they would 
satisfy the identity 

erf(efc/v^) +erf(e^_,+i/v^) = 1, 

which is clearly inconsistent with the exact relation (14.41) . However, the 
slightly modified formula 

a - V2 erf- (^^^^) (4.13) 

differs from (14. 120 by a term of order 1/A^ and is fully consistent with the 
relation (14.41) . Although both (14.121) and (14.131) provide an excellent approx- 
imation to the chain sites for large N, the latter equation is always more 
accurate than the former, and can be used to estimate C,k with remarkable 
precision even for relatively low values of A^, cf. Fig. 14.31 □ 

In the rest of this section we shall analyze several relevant features of 
the site distribution using the above heuristic approach. It is certainly of 
interest to determine whether the position of the last spin tends to infinity 
as — >■ oo, since according to our interpretation of the chain's geometry the 
number I^n/t^ is the radius of the circle on which the spins lie. According to 
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k 



Figure 4.3: Sites coordinates and their continuous approxima- 
tion flTOj) for = 20 spins. 



Eq. fl4.13p . for large the last spin's coordinate is approximately given 
by 

~ erf-i (l - ^) , (4.14) 

so that should diverge as — oo. Of course, this assertion should be 
taken with some caution, since in Eq. f l4.13p the argument of the inverse 
error function is correct only up to terms of order 1/A^. In order to check the 
correctness of the approximation (14.141) . we use the asymptotic expansion of 
erf~^(M) for u — > 1 in Ref. [26] to replace (14.141) by the simpler formula 

^ V2r/-logr/, (4.15) 

where 

^ = (^) • 

As can be seen in Fig. 14.41 the approximate formula (14.151) qualitatively 
reproduces the growth of when A^ ranges from 100 to 250. A greater ac- 
curacy can be achieved by introducing an adjustable parameter in Eq. (I4.14p 
through the replacement 1/A^ — ^ a/N, so that in Eq. (14.151) r] becomes 

/ A^ \ 

r/ = log^ . (4.16) 
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In Fig. 14.41 we have also plotted the law fl4.15p -( H:.16l) with the optimal 
value a = 0.94, which is in excellent agreement with the numerical values of 
for iV = 100, 105, . . . , 250. 

The last property of the spin chain (14.11) that we shall analyze in this 
section is the dependence of the coupling between neighboring spins on their 
mean coordinate. Calling 

we shall now see that when > 100 the Gaussian law 

Ar2 - 2 

e-«'= (4.17) 

holds with remarkable precision, cf. Fig. 14. 5[ Indeed, if x = x{k) denotes the 
RHS of Eq. dUSD we have 

2k = Ned (^-^) + + 1, 



so that 

dx \/27r 
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G 




Figure 4.5: Plot of the coupling between the spins k and k + 1 
versus their mean position, fitted by the Gaussian fl4.17p . for = 
100 spins. 



is of order 1/A^. Hence, up to terms of order 1/A^ we have 

^2 



Since (up to terms of order 1/A^^ 



\dk 



27r 



x{k) ~ - 



x[k--)+x[k + - 



(4.18) 



Eq. (14.171) follows from (I4.18P replacing k hj k + ^. 

A brief comparison of the previous properties with those of the PF chain 
(I1.28P is in order. For large A^, the density of sites of the PF chain (normalized 
to unity), that is the density of zeros of the A^-th Hermite polynomial, is 
asymptotically given by [36] 



pN[X) 



7CN 



V2N 



x^ 



(4.19) 



The last site (n of the PF chain grows with much faster than the corre- 
sponding site of the chain (14. ip . for the largest zero of the A^-th Hermite 
polynomial behaves as ^/2N + 0{N~^^^); see, for instance, the recent pa- 
per [56]. 
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4.3 Quasi-exact solvability of the NN chain 

In this section we shall explore the solvability properties of the spin chain fl4.ip 
by exploiting its connection with the QES model Hq defined in f l3.2ap . As in 
Section we shall rescale the constant u in the latter Schrodinger operator 
and its scalar counterpart H^^ and remove the ground state energy, defining 

H = {Hq — Eq) \^=a , H'''' = {Hq" — Eq) \^=a ■ 

In terms of the matrix multiplication operator 

h(x) = ^(xi - - Si^i+i) , 

i 

the connection between H and H^" can be simply written as 

H = H''' + 2ah , (4.20) 
whereas the spin chain Hamiltonian (14.11) is recovered as 

H = h(0. (4.21) 

Recall that the ground state function of the scalar Hamiltonian H^" can 
be written as /i(x) = e"^*^''\ cf. Table ISTTl and Eq. (14.31) . We shall also denote 
by 

<^o = ||/i|rVeLi(Co)nC°°(Co) 

its normalized square. The following result is easily established by repeating 
the proof of Lemma 12.221 

Lemma 4.5. ipo converges to the delta distribution supported at ^ and decays 
exponentially fast away from ^. 

One would be tempted to beheve that, mutatis mutandis, the discussion 
of Section (12. 5p would provide an appropriate bridge to link the algebraic 
states of H and H. However, the fact that the dynamical models H and 
H^'^ are only quasi-exactly solvable, and thus only a proper subset of their 
spectra is known, makes it impossible to compute their partition functions 
from their algebraic energies. Hence Lemma 12.261 which constitutes the core 
of the application of CS models to the study of spin chains, cannot be invoked 
in this case. 

In the rest of this chapter we shall develop a finer version of the freezing 
trick that will enable us to circumvent this limitation. We shall obtain a 
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number of eigenstates of the chain (which we call algebraic), valid for any 
choice of and M. As mentioned in the introduction, these algebraic states 
do not span the whole spin space, the chain Hamiltonian fl4.1l) thus inheriting 
the QES structure of the dynamical model. 

In the following lemma we prove the key convergence result ("freezing 
trick") that shall enable us to compute the algebraic eigenvectors of H in 
closed form. This result is most easily stated in terms of the function /t = 
that is, the normalized ground state of H^'^. Moreover, given a 
scalar function ip G L^(Co) and a spin state $ G L^(Co) ® S, we shall use the 
notation 

((^,$) = / ^$ 
J Co 

to denote their S- valued scalar product. 

Lemma 4.6. Let \E' be an algebraic eigenfunction of H with energy E (cf. 
Theorem 3.11^ . Assume that there exists an eigenfunction ip of H^^ and a 
Tj-valued polynomial F G C[x, a~^] ® S such that 

/i-^(^ -^F) = C(a"^) , (4.22) 

and suppose, moreover, that these eigenfunctions can be normalized so that 

/i-^^ G C[x,a-^] ® S, /i^V e C[x,a-^] . (4.23) 

Then the limits 

c = hm , X = lim F{^) , E = lim ^ 

a— >oo UlcJ a— »oo a— >oo ZCt 

exist, and 

c(H-E)x = 0. (4.24) 

Proof. The existence of the above limits trivially follows from the polynomial 
dependence of fiip and F on and the fact that E is linear in a. Further- 
more, the self-adjointness of H'^^, Eqs. f l4.21l) - P:.22l) and Lemma (4.51 readily 
imply that 

(/i, H'''^) = {H'^'fi, ^) = , 

(/i,h^)= / ipo[ix-^^hF + Oia~')]=cHx + Oia-^). 

Jco 

If we now use Eq. f l4.20p and Lemma 14.51 to write 

(/i, W'^) + 2a(/i, hv]/) = (/i, H^) = Fax + 0{1) , 

substitute the former equations into this identity and take the limit a — > oo, 
we readily derive Eq. fl4.24p . □ 
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Corollary 4.7. If cx 7^ 0, then x is an eigenvector of H with energy E. 
Remark 4.8. Since, by Lemma [4.51 

c = lim (/i, Tp) , 

a^oo 

the previous corollary is of no practical use unless ip can be taken as the 
normalized ground state /i, in which case c = 1. 

We are now ready to prove the main result of this chapter. 
Theorem 4.9. // |s) G S and |s') G S', the states 



Xo 


= Xo(|s)) = 


A|s), 




(4.25a) 


Xi 


= Xi(|s)) = 


i 




(4.25b) 


X2 


= X2(|s')) = 


i 


-i)5^6e,is:+), 

i<j 


Si2|s') = |s') , (4.25c) 


X-i 


= X3(|S)) = 




5i2|s) = -|s), (4.25d) 



satisfy the equations 

(H-2)x. = 0, ^ = 0,1,2,3. 

Proof. The proof consists in the application of Lemma 14.61 to appropriate 
triples (\l/,?/',F), where in fact ip = fi (and hence c = 1) in view of the 
previous remark. 

In the first place, it is easy to show that xo ^ ker H directly. However, 
we prefer to prove it using the same freezing trick argument as in the other 
cases. To this end, consider the functions 

^(0) 

^ = ^ = /iA|s), F = A|s), 
11/^11 

which trivially satisfy the hypotheses in Lemma [4.61 with E = and ip = fi. 
Since in this case 

E = 0, X = A|s), 

from Corollary 14.71 it follows that X05 if non- vanishing, is an eigenvector of H 
with eigenvalue 0. 
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The next case is slightly more complicated. Let us set 
11/^11 

It is not difficult to show that (\E', fi, F) satisfies the hypotheses in Lemma 
with E = 2a. Since E = 1 in this case, Corollary 14.71 shows that the state 

X= lim F(|) = ^e.|s.) 

i 

is either zero or an eigenvector of H with energy 1. Here we have used the 
identity Ylii^i — 0' "^f- Corollary 14. 5[ 

The other two states are obtained in a similar manner starting with ip = fi 
and the functions 



^ ^ ^02 +1 ijy^02 ^ ^^^(2) ^ _ ^)^(2) _ 2iVx$(i) + {N + 2)x2$W) 
I A'' 



F = Z, E = Aa, |s')gS', ^lajs') = |s') 



and 

$(3) _ 4^(1) _ 8^(0) 

^03 ^^03 3 03 



$(3) _ 2x<l>(2) + ^$(1) + X - ^ ] $ W 



11/^11 

F = -, E = 6a, ^i2|s) = -|s) 
/i 

in each case. Indeed, with the former set of functions we immediately arrive 

at 

i i<j 

while in the latter case we obtain 

i<j i 



Using Proposition 14.21 the result follows. □ 

Remark 4.10. All the algebraic energies of H being integral multiples of a, it 
follows that any eigenvalue of the chain constructed from them via a freezing 
argument must be an integer. 



CHAPTER 4. NN SPIN CHAINS OF HS TYPE 



107 



Remark 4.11. The linear combinations of eigenfunctions with energy 2 and 
3 used in Theorem 14. 9 1 are not casuaL Should one start, e.g., with 

a\\n\\ 

for which F = 0{1) and c = 1, one would be lead to 

2N 



X 



N 



$(0) 



The above linear combinations have been chosen so as to suppress the dom- 
and §g) 



inant parts in \1>q2'' and "^m, which respectively are proportional to x'^ and 



It is also clear that the state X2 (resp. Xs) is zero whenever the spin 
vector |s) (resp. |s')) is not symmetric (resp. antisymmetric) under 5*12 , which 
accounts for the condition on |s) and |s') imposed in Eqs. (]4.25bll and ( I4.25dl) . 

It is natural to wonder whether additional eigenvectors of the chain can 
be computed using some different linear combinations of the algebraic eigen- 
functions of H listed in Theorem 13. 17[ Nevertheless, since 

cx = lim F(^) = hm — — , 

it is clear that an algebraic eigenfunction in the linear span of "^[{l and 
can only give rise to a nonvanishing vector cx when / = and m = j. 
Indeed, = a'^ jlLf^^\ar^) e ;uC[x,a-i], with /3{0) = 1 - N{a + i), is an 
eigenfunction of H^'^ with energy 4a/, which implies that 

lim a-^L;^~\aN^) = V/ > 



a-s IICIP = N, proving the former statement. The latter condition follows 
directly from the factor x™"-' present at the above eigenfunctions of H. This 
discussion strongly suggests (albeit does not rigorously prove) that one can- 
not obtain any additional eigenvectors of H using the technique developed in 
Lemma 14.61 and Theorem 14.91 



4.4 Qualitative properties of the algebraic spec- 
trum 



In this section we shall analyze which portion of the spectrum of (14.1 p is 
algebraic and how the algebraic eigenvalues are embedded into the whole 
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spectrum. Regarding the first point, it is convenient to start with the follow- 
ing states $i e L^(Co) ® S, whose evaluation at ^ yields the chain eigenvec- 
tors Km -. 

<|.o(|s)) = <|.W(|s)), 
<|.i(|s)) = <|.«(|s)), 
<|.2(|s)) = <|.(2)(|s)) + (iV-l)$(2)(|s)), 

$3(18)) = <|.(=')(|S)) - 2X^^'\S) + — $«(|S)) . 

It may be easily shown that if |s) is symmetric then $i(|s)) and $2(18)) 
vanish at x = ^. Likewise, if |s) ~ then $2(|s)) is also zero when x = ^. 
For this reason, we shall consider the following spaces: 



Vo = A(S) , Vi = ($i(|s)) : |s) e S/A(S) 
V2 = ($2(|s)):|s)g(s7A(S))/~), 
V3 = ($3(|s)):|s)gS, Si2|s) = -|s)\. 



Proposition 4.12. The dimensions of the spaces Vi are given by 

dimVo=(^ ^ j, (4.26) 

/N + 2M - 1\ 

dim Vi = dim V2 = (AT - 1) ( ^ ^ j , (4.27) 

fN -l\fN + 2M -2\ , , 

dimV3=( 2 )( ^ )■ (4.28) 

Proof. First of all, the dimension of Vo is simply the number of permutations 
with repetitions of A^ elements from 2M + 1. 

Consider next the space Vi. If two basic states |s) and |s') differ by a 
permutation of the last A^ — 1 spins, then $i(|s)) = $i(|s')). Hence the 
dimension of Vi is given by 



(2M + 1) 



AT + 2M - 1\ /N + 2M 
2M ) ~ [ 2M 



where the first quantity is the number of possible choices for si, the second 
one is that of symmetric states of A^ — 1 particles, and the last one, that of 
symmetric states of A^ particles. 
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The dimension of V2 coincides with that of the space (S'/A(S))/~, a 
basis of which was computed in Propositions 13.101 and 13.131 Using the pre- 
vious notation we can say that, for a given spin content C 
{-M, ...,M} with muhiphcities (z/i, . . . , G N" (with i^j = A^), 

there are n — 1 independent states in (S'/A(S)) /~ . Hence 



min{2M+l,Ar} 



"""^"^"^"'^'^ /2/Vf_|_l\ 

dimV2= - 1) card(<P;v n M") ( J, 

n=l \ ^ / 

where the binomial term accounts for the different choices of {s^, . . . , s"}, 
and the cardinal for all the possible values of the multiplicities. Noting 
card(<P7v n N") = {^zD and using the identity 



n\/ m \ I m + n 
k)\j + k) \m-j 



E 

A;=0 

we arrive at the desired formula. 

It only remains to prove Eq. (14.281) . To this end, let us concentrate on 
the spin vectors Sq with a fixed spin content {s^, . . . , s"}, (z^i, . . . , z/„). It is 
clear that V3(So) = ($3(|s)) : |s) G Sq) is spanned by the (2) states 



V s , s 



Z/,Z/,$3(|s^S^---)-|s^S^---)) 



We shall shortly show that, for all i, 



v{s\s^) = 0. (4.29) 



By antisymmetry, the double sum ^ v(s*, s-') vanishes, so there are at 
most n — 1 independent relations of the above form. It can be verified that 
{'v{s\s^) : 1 ^ i < j < n} is in fact a basis of V3(So), so the dimension of 
V3 is obtained from (ii) by substituting n — 1 for ("2 ^) . 

Finally, let us prove (I4.29p . A first observation is that 

r'$(i)(|sV---)) - Ara;$(')(|sV---)) = -J^'^k^^^'Hls's'' ■ ■ ■)) . 



k 



Hence 

A^ 
2iy, 



^^v(.\.^-) = Ar5^z.,$(3)(|.V---)) 



J2 i^ji^k^^''\\s's'' ■■■)) + 5^ z/,z/,,$(=')(|s^s'^ ■■■)). (4.30) 
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The first two terms of tfie RHS obviously cancel, and tlie last one vanishes 



This proposition, which is interesting in its own right because of its con- 
nection with the dynamical model fl3.2al) . has the following consequences for 
the spin chain (14.11) : 

Corollary 4.13. For each i = 0,1, 2, 3, the number of algebraic eigenvectors 
with energy i is bounded by dimVj. 

Corollary 4.14. The ratio 



total number of states 

tends to zero as N oo. The whole kernel of H is composed of algebraic 
eigenvectors. 

Proof. The first part of the statement follows by direct comparison with the 
total number of states (2M+ 1)^. The second part is due to the fact that H 
is a linear combination with positive coefficients of the nonnegative operators 
1 — which implies that ker H = A(S). □ 

These observations are corroborated by our numerical simulations: the 
number of algebraic levels increase as M and become larger, but never- 
theless they span a lesser and lesser part of the spin space. The kernel of the 
chain, e.g., shows clearly this behavior. But most remarkable is the follow- 
ing conjecture regarding the embedding of the algebraic eigenvalues in the 
spectrum. 

Conjecture 4.15. The algebraic eigenvectors fl4.25ap - fl4.25cp span the three 
lowest levels of the spin chain fl4.ip . When M ^ 1 the fourth algebraic energy, 
E = 3, is however the fifth lowest level of the chain, and its whole eigenspace 
is spanned by the states of the form f l4.25d|) . The algebraic energies are 
singled out in the spectrum of H as being the only integer ones, i.e.. 



We have numerically diagonalized the matrix representing the Hamilto- 
nian H for M ^ 2 and up to 12 particles, and our results fully support 
the above conjecture. Another important observation is that the standard 
freezing trick relation between the energies of the spin chain H and of the 
dynamical models H and H^'^, namely 



by antisymmetry. 



□ 



number of algebraic states 



spec(H) nZ = {0,1,2 



3}. 



E 



lim 



E-E'' 



'SC 



a- 



■oo 



2a 
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and the fact that H has noninteger eigenvalues imply that the spin dynamical 
model fl3.2ap has noninteger energies, and hence that some of its eigenfunc- 
tions are not among those listed in Theorem I3.17[ 
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Chapter 5 

Conclusions and perspectives 



In this thesis we have studied three QES models of CS type presenting short- 
range interactions and developed a rigorous treatment of the CS/HS corre- 
spondence, valid for both ES and QES Hamiltonians. In a nutshell, the main 
conclusions and results that one should extract from this dissertation are the 
following: 

(i) On the BCn HS spin chain: We have computed the spectrum of the 
trigonometric BCn Sutherland model by expressing its Hamiltonian as 
the sum of squares of an appropriate set of commuting Dunkl operators. 
The BCjsi HS chain is obtained from this system by "freezing" the 
particles at their classical equilibrium positions. We have proved a 
rigorous version of the freezing trick which has enabled us to compute 
the partition function of the chain in terms of those of the corresponding 
scalar and spin models. We have resorted to numerical computations to 
explore the level distribution of the chain, finding out that the energies 
accurately follow the Gaussian law. A constant magnetic field was 
added to the picture in a rather straightforward way. 

(ii) On NN spin models: We defined three QES spin models with near- 
neighbors interactions that are related to NN versions of the sum of 
the squares of appropriate sets oi Dunkl operators. Furthermore, 
we have obtained infinite fiags of finite-dimensional spaces invariant 
under these models, and computed several families of eigenf unctions 
in closed form by diagonalizing the restriction of the Hamiltonians to 
these spaces. All the algebraic solutions, whose energies are integers, 
can be expressed as the ground state times a symmetric function. We 
conjecture that some of their eigenfunctions are not of this form, and 
that their corresponding eigenvalues are not natural numbers. 
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(iii) On NN spin chains: We have defined a spin cliain presenting nearest- 
neighbors interactions that is related to one of the previous NN models 
by a freezing trick argument. We have discussed its site distribution 
in detail, both analytically and numerically, and obtained a number of 
algebraic eigenvectors and energies using a refinement of the techniques 
developed in Chapter [2] that also applies to QES systems. We have 
raised the conjecture that the only integer energies of this chain are 
0, 1, 2, and 3, and that the eigenstates corresponding to these energies 
are all algebraic. 

The conjectures scattered throughout the text do not exhaust all the open 
questions that we find worth exploring in connection with our work. In what 
follows, we shall briefly mention some wide research lines that in our opinion 
would be of interest, and that we intend to follow in the near future: 

(i) The structure of the (non- algebraic) eigenfunctions of Hq (and other 
QES operators): From the previous discussion, it looks unlikely that 
all the eigenfunctions of Hq admit the ground state as a Jastrow factor. 
Nevertheless, this does happen in the theory of CS models, and the 
functions that are explicitly constructed in the usual QES models are 
generally of this form. A thorough study of this phenomenon is, to 
the best of our knowledge, still lacking, and would probably shed some 
light on the connections between Schrodinger operators and orthogonal 
polynomials. The particular case of Hq would be of special interest. 

(ii) HS/CS and symmetric spaces: A powerful procedure for constructing 
CS models with spin starting out with a symmetric space of nonpos- 
itive curvature has been recently developed by Feher and Pusztai, as 
described in Section 11.51 It should not be difficult to extend the tech- 
niques described in this thesis to solve the associated spin chains. 

(iii) Statistical properties of the spectrum of CS models: The remarkable 
statistical properties of the BCn HS chain occupied a significant por- 
tion of Section 12.61 It is difficult to obtain pen and paper results on 
these aspects, but the rigorous literature studying these topics is al- 
ways increasing (see [38,49] and references therein). Further study of 
the topics mentioned in this section and extensions to the case of NN 
models would make a valuable addition to the existing literature. 

(iv) Integrability vs. solvability in spin chains of HS type: A satisfactory 
analysis of the integrability of the chains of HS type has only been 
accomplished in a few instances. A natural approach to this problem 
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would be to define a set of integrals of the chain using the freezing 
trick, but in general it is not easy to grant the linear independence of 
these conserved operators. If one managed to overcome these technical 
problems, this method (which has already been exploited in the study of 
some particular models) would provide a clean proof of the integrability 
of the chains associated with integrable spin models. 
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